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Reference Abstract 

Documentation of three model versions: the Douglas-fir tussock 
th population-branch model on (1) daily temporal resolution, (2) insi 
mporal resolution, and (3) the Douglas-fir tussock moth stand-outbre< 
del; the hierarchical framework and the conceptual paradigm used are 
scribed. The coupling of the model with a normal-stand model is dis- 
ssed. The modeling convention and algorithm used are also discussed, 
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INTRODUCTION 

This paper documents the mathematical form and computatioi 
of the Douglas-fir tussock moth (DFTM) stand- outbreak model, 
is constructed to simulate outbreaks of Douglas-fir tussock mo 
pseudotsugata (McDunnough) ) , a defoliator of Douglas-fir (Pseut 
menziesit var. glauoa (Beissn.) Franco), and the true firs (Ab* 
in western North America. The details of the model's backgroui 
historical development, and the basis for formulation of its p: 
structure and functional forms will be covered in another pape: 
documentation of model behavior. 

The conceptual structure of the DFTM model provides a per! 
on the relation of the various models developed. Three versioi 
DFTM outbreak-population model were developed. Version 1.2, pi 
here, is at daily resolution and differs little from the origii 
(VI. 1) conceived in 1974. The second version presented (V2.2) 
originally developed (V2.1) to reduce computing costs and necei 
only minor alterations. These two versions follow the insect < 
branch through 1 year. This branch model became the operation; 
the stand-outbreak model (V3.1). Details of the coupling relai 
between the normal-stand model and the DFTM stand-outbreak mod* 
the overview of the general model structure. 

A summary of the FLEX processing algorithm and notational 
precedes detailed documentation of the models. Adjustments to 
current algorithm were necessitated by the DFTM stand-outbreak 
structure, and these are detailed next. The computational frai 
which the models are embedded is discussed in these two sectioi 

A short discussion, preceding the details of each model v< 
gives the relations among the three model versions and how the: 
from each other. The detailed FLEXFORMS, summarizing model st: 
and parameterization, follow. I/ Next are the specifics of fill 
used in the third model version. This includes a more detaile< 
of the computational order discussed in general in another sec 

Last is a description of the stand-postoutbreak model. T" 
the filtering mechanism that assigns mortality and growth impa 
stand of trees. These result from defoliation of the model br. 
transformed by the whole-tree-defoliation equation. 



A hierarchical structure is conceptualized for the general model: 
ich level is a dynamic unit of the level above (figures 1 and 2). 
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Figure l.--How various models 
assemble into the hierarchical 
model. A region (R) consists 
of a number (L) of forests (F), 
a forest consists of a number 
(M) of stands (S), and each 
stand consists of a number (N) 
of trees (T) . 
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Figure 2 . --Hierarchical or- 
ganization of the DFTM stand- 
outbreak model relative to the 
normal-stand model. 



A. Regional model: The socioeconomic context of the forest. 
Resolution: 1 century; not implemented. 

B. Forest model: An explicit assemblage of ecotype-stand units. 
Resolution: 1 rotation; not implemented. 

C. Stand model: 

(a) Normal-stand model: An assemblage of trees of varying 
size and species, in a specific ecological and management 
regime. 

Resolution: 1 decade; not implemented. 

(b) Stand-outbreak model: An assemblage of tree classes 
prescribed by the effective resolution of knowledge of stand 
conditions and insect population levels. The conceptual 
structure of an outbreak is patterned after the observed 
sequence of population "phases" of Wickman et al. (1973) . Tt 
successive years of an outbreak are parameterized according t 
phase-specific mortality vectors and other population propert 
Resolution: Duration of outbreak, 4 years. 



defoliation into mortality, top-kill, and reduction of 
for each tree class and covers delayed mortality and r 
management practices. 

Resolution: Duration of postoutbreak period, 6 years, 
(d) Stand-annual model: The potential of integration 
tree classes at the end of each year effectively impos 
stand-annual structure on the assemblage of tree class 
Redistribution of eggs over the tree classes in the st 
determined in this model. 
Resolution: 1 year. 

D. Tree model: 

(a) Normal- tree model: The rules of normal growth an 
survival of a tree, of specific species, size, and age 
of specific ecological and stand regimes. 
Resolution: 1 year; not implemented. 

(b) Tree-class model: A class of trees is represente 
typical tree. Defoliation of the typical tree over a 
is derived from the branch model. 

Resolution: 1 year. 

(c) Tree-outbreak model: Normal growth increments an 
survival are modified by the degree of defoliation as 
mined from the branch model. The tree-outbreak model 
simply a branch model, coupled or linked over seasons, 
outputs translated at the end of the outbreak into the 
of defoliation on tree mortality, top-kill, and growth 
model does not explicitly appear, being subsumed by th 
outbreak model. Effects of defoliation are relegated 
stand-postoutbreak model. 

Resolution: Duration of outbreak and postoutbreak per 

E. Branch model: Contains the rules of growth and survi\? 
tussock moth, branch defoliation, and the production c 
foliage and loss of old foliage by the branch; specif i 
ratio of eggs to foliage, for nominal ratio of new to 
foliage, and for species of tree and year (phase) of c 
Standardization of the total biomass of foliage on the 
branch (100 grams for Douglas-fir and 144 grams for gi 
(Abies grandis (Dougl.) Lindl.) allows specification c 
two variables, percent new foliage and eggs (per total 
of foliage). 

The mechanisms of the branch model are conceptual 
a daily resolution, with runs spanning one season. Ir 
conditions for a season run are identifiable from the 
nating condition of a season run which provides the c 
to couple several seasons (years) together. In the dc 
resolution branch model, days are occasions. 



The postoutbreak period has been set at 6 years so that the 01 
postoutbreak periods will sum to the resolution of the normal-stand mi 



eixiciency in computing, tne aaiiy resolution moaei 
was collapsed into an instar resolution model. This changes 
little and reduces the number of occasions in a season. 

We can now identify several levels of temporal resolution: outbreak 
aining several phase-specific years) , year (containing several 
ions), and occasion. 



IDENTIFICATION OF COUPLING VARIABLES 

In this outline of the couplings between the different models, 

c numerals correspond to those in the numbered diamonds in figure 2. 

3utputs of the normal-stand model and inputs to the classification 
nodel: Those state variables and parameters of the normal-stand 
nodel needed to assign trees of the stand to the classes of the 
stand-outbreak model. 

Dutputs of the classification model and inputs to the stand-outbreak 

nodel: 

L. The number of tree classes in the stand; 

ii. For each tree class, represented by a "typical tree": 

a. The number of trees in the class; 

b. The foliage composition of the model branch: Biomass of 
foliage and percent new foliage; 

c. Egg density; 

d. Species; 

e. The weighting factor to account for relative contribution 
of a tree in the class to stand foliage biomass and 
insect population density; 

Lii. Ecological or site-specific model parameters. 

Outputs of the stand-outbreak model and inputs to the stand-postoutbr 

nodel for each class of trees in the stand: 

L. Defoliation of the model branch after each season during the 

outbreak; 
Li. Tree defoliation at the end of the outbreak. 

Dutputs of the stand-postoutbreak model and inputs to the 

treclassif ication-ef fects model, for each tree class: 

L. Direct mortality; proportion of trees that die as a direct 

result of defoliation; 
Li. Secondary mortality; proportion of trees that suffer secondary 

mortality in the postoutbreak period, identified as: 

a. Bark beetle mortality; 

b. Other mortality caused by the outbreak; 

c. Noncausal mortality; 

Lii. Top-kill; proportion of trees in each top-kill subclass. Each 

subclass is defined in terms of height reduction. An additional 
output, the proportion of the nominal diameter- growth increment 
actually achieved is also given; 
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a. Proportion of the nominal height-growth increment actiu 
achieved ; 

b. Proportion of the nominal diameter-growth increment 
actually achieved; 

v. Factors for delayed reduction of growth beyond the postoutbi 
period. 

Outputs of the classification model and inputs to the reclassi- 
f ication-ef fects model: The record of which trees were assigned 
the various classes. 

Outputs of the reclassif ication-ef fects model and inputs to the 

normal-stand model: The allocation of mortality factors to the 

trees of the stand and the updating of growth increments to covei 
the outbreak period. 



ABSTRACT OF THE FLEX CONVENTION AND ALGORITH 

The conceptual basis of FLEX is the identification of several ki] 
f variables and functions in the context of a general algorithm for 
ynamic behavior. Identify: 

x, the vector of state variables; 

y, the vector of output variables (functions of state var: 

z, the vector of input variables; 

b, r, two vectors of constants (-Rodel parameters). 

he update algorithm is a simple difference equation 

= x(k) + A(k) 



ith A(k) constructed as a sum of internal variables that are identi: 
s f function values. These functions account for fluxes (flows) bet) 
tate variables. Flux from x. to x. is f . . ; 

f = f (k,z,x,g,s,b,r) ; 

-LJ J.J ~ - ~ ~ ~ ~ 

nere g and s are vectors of functions; f-^ is a source or sink for 2 
spending on the sign of fj^. The g functions are intermediate funct: 
iternal variables that are used extensively, 

gi = g (k,z,x,b,r). 

le s functions, or special functions, are included for further model: 
id programing flexibility; these are implemented as FORTRAN-?/ subroul 
ad are therefore much more general. Details of how this process is 
nplemented are given in the FLEX users 1 manual (White and Overton 
377). 

FORmula TRANslating system- -a language used to express computer pro 
f arithmetic or logical formula. 



and from this point forward uses the values of the variables at time ] 
to compute the values of the state variables at time k+1. The procesi 
uses the values for time k to compute the outputs for time k+1. This 
accounts for actions taking place during occasion k+1. 

| -< ....... occasion k > | <- - occasion k+1 - > \ 

time time time 

k-1 k k+1 

In, a terminal (or FLEX) module the algorithm may be summarized a; 
follows: 

(1) Receive (or calculate) inputs and store the values of z 
using the z functions. 

(2) Calculate and store the values of the g functions. 

(3) Calculate and store the values of the f functions. 

(4) Construct the update vector, A. 

(5) Update the state variables. 

(6) Calculate and store (generate outputs) and values of y usin; 
the y functions. 



k time step k may be subdivided into q sub intervals indexed by k T = 
L, ..., q - 1. Then define: 

x(k) = x(k, 0), 

x(k, k T + 1) = x(k, k') + A(k, k' ) , 

x(k + 1) = x(k, q); 

f is added to the g and f function arguments. This capability impli< 
an additional step, 5a, in the above sequence: 

(5a) Increment k T and compare to q where q is the resolution 

ratio. If less than q, return to step 2, otherwise proceed 
to step 6. 

Processing of a ghost (REFLEX) module, one that has control over 
submodules, requires a modified procedure; steps 3, 4, and 5 are repL 
Dy step 3a: 

(3a) Call and process subsystems (submodules) , which in turn upd; 
the values of the state variables in ghost. Then continue 
with step 5a. 

The current FLEX algorithm assumes that a division of one time-s- 
into q subintervals within a terminal (FLEX) module is being done to 
implement a solution to a set of differential equations, and according 
:>nly 1/q of the update (A) value is added to each variable. In modul< 
$1 of version 3, this capacity is used differently; an event structur< 
is assumed and a full update on each variable is desired. Because 
q = 10 (representing the 10 occasions in a season) , each f function ij 
nultiplied by 10 to nullify the automatic multiplication by 1/q (Whit< 
and Overton 1977, p. 21 and 32). 



computational convenience internal to the FORTRAN code of a function 
(e.g., a g function). These are available to no other functions and 
tfill not appear in the list of explicit state variables. 

VECTOR PROCESSING IN FLEX2 

The current version of the DFTM stand- outbreak model makes use o 
some of the conceptual structure of the FLEX paradigm that is not 
explicitly a part of the current algorithm used by the model processo 
is detailed in the FLEX users manual. The following explanation will 
clarify a newly adopted notational convention, its conceptual meaning 
the paradigm, its place in the algorithm, and its implementation usin 
the current version of the processor. 

In the FLEX paradigm, x, the vector of state variables, may cont 
elements that are themselves vectors usually when a class II system 
)eing considered, a system in which each subsystem has the same struc 
md function, as in a population of individuals that have the same ru 
3f behavior. Thus, Xj_ may be an element of x. For example, x^ might 
represent the biomasses of each of the individuals in a population, w 
in element of Xj[ the biomass of a particular individual. The notatio 
convention is established that x is the vector of state variables; x-^ 
i vector element of x; and x^ is a scalar element of x^. 



Only the homogeneous case is allowed, so that each member subsys 
Df the class 11 system has the same state variable structure. A g fu 
:hat operates on x^ is a vector valued function. g is the vector of 
Intermediate functions with gj_ an element of this vector. g^ is a sc 
/alue of g-j^ where the caret notation is used to signify an actual val 
Df the function. When g^ operates on the vector x, it may produce a 
/ector g.j_. The values generated by operation of g on x form a matrix 
corresponding to the ith element of g and the jth~member of the popul 
Chen the notation g.^ denotes the column vector of values for all g f 
cor the jth population member, and the notation g^. represents the ro 
/ector of values of the ith g function for all population members. A 
similar notation is extended to the function vectors z and y; b and r 
parameter vectors may also have vector elements. 

Introduction of the ability to cycle through a single subsystem 
specified number of times, with each cycle processing in turn a membe 
Df a population of identically structured subsystems, required modifi 
cation in the FLEX algorithm. Instead of incrementing k' automatical 
ifter each pass through the subsystem, counter k !f is incremented afte 
*ach pass and then checked against p. The population number, p, is s 
)efore a simulation run. If the subsystem has not been processed p 
:imes, then the g functions are called and the subsystem is processed 
again. The values g ^ are calculated, the updated values of x^j are 
md k" and j are incremented. After cycling p times, the vector elem 
Df g and x have been updated; k' is then incremented and the next pas 
through the g functions will fix g.., before the first call to the sub 
Eor that time step. 



ite variables. A, is updated, by column vectors and tne values u are 
.culated by column vectors. It implies, however, that internal storag 
nations for the values g^^ , x-j^ , etc., are within the processor. Thes 
nations do not exist in the current processor version. Vector values 
Jt be stored on external files that are accessed, manipulated, and 
lated by the various functions and subroutines that use them. 

The following procedure is used to implement this algorithm feature 
the current processor version. At k' =0, the g functions are called 
usual. In a ghost module, each g-j_, may operate in a different mode. 
5^ may calculate g^. at this time and store it on an external file. P. 
T calculate only a scalar value and have it available in the standard 
7 for use by other gj/s. Or a g^ may be used to process external file 
Ltialize the state variable and input variable matrices, and store 
im on external files for subsystem access. Each time the subsystem is 
)cessed, the functions within the subsystem reads ghost's g 4 from the 
_e. The g functions in the subsystem generate only that subsystem's 
_umn vector g.-j, and these values are stored in the processor and are 
:essed in the standard way. The subsystem's updated values of x . are 
>red on a different external file, not in the ghost module. On comple 
subsystem processing, the functions in the ghost module may then 
:ess those files as necessary. 

The above procedure introduces a new level of complexity in FLEX 
>graming. Much of this complexity could be avoided by introducing a 
7 level in the hierarchy, but this entails difficulties as well as 
ireasing processing time and cost. The approach used here is relative 
expensive and easily implemented. 



CORRESPONDENCE OF THE THREE MODEL VERSIONS 

The model for interaction of the tussock moth larvae and host 
.iage was originally conceptualized and developed in terms of daily 
;olution mechanisms. Table 1 defines the daily resolution occasion 
:ucture of version 1.2. This model was restructured into the instar 
;olution of table 2 for version 2.2. Only one change in the conceptua 
ructure is necessitated by this change in resolution. The model 
)wth and feeding of the tussock moth is altered in the fourth instar 
accommodate the convolution of larval growth with population survival 
ir an instar. This is accomplished by assuming that the growth rate 
constant throughout this instar, at the average of the two growth 
;es used during this instar in version 1.2. This changes the accumula 
:val biomass, but only slightly and only during the fourth instar. 

The organization of the models by subsystem and function identifica 
shown in table 3. The subsystems are those of the branch model at 
:h daily and instar resolution. The g functions, although representin 
: same processes, do change form from version 1.2 to version 2.2. The 
:ails of the functional forms are given in the FLEXFOKMS that follow. 



del of module S-^ of version 3.1 (see figure 3). The only changes 
cessary are: 

(i) some "normal tree 11 parameters are being brought in as z inpul 
ther than being internal b parameters; 

(ii) the time index is changed from k to k ? ; 

(iii) the overwintering activities of occasion 10 of table 2 take 
ace in the ghost module SQ of version 3.1. The details of these 
anges are presented in the version 3.1 FLEXFORMS. 

Table 1 Douglas-fir tussock moth population-branch model, daily 
resolution form, Version 1.2. Occasion/time structure: 
note that occasion t is the interval before time k; 
k = t-1 during occasion t 

Occasion Description 

Initiation Bud burst and shoot elongation complete. 

1 Egg hatch and establishment of first instars. 

2 First day of feeding, growth, and mortality. 



10 Ninth day of feeding, growth, and mortality. 

11 Last day larvae are in the first instar: feeding, growth, 
and mortality. 



21 Last day larvae are in the second instar: feeding, growth, 
and mortality. 



31 etc. 



41 etc. 



51 Last day larvae are in the fifth instar: feeding, growth, 
mortality, and male pupation. 



61 Last day of sixth instar: feeding, growth, mortality, 
and female pupation. 

62 Pupal mortality and adult female emergence. 

63 Eggs laid, aging of new foliage, and new foliage 
potential set. . 
Redistribution of insects over the stand. i/ 

64 Overwinter mortality of foliage. 
Overwinter mortality of eggs./ 

This is not operative in a single-tree run, only in the stand- 
outbreak model. 

Overwinter mortality accounts for all the mortality factors to 
which the egg masses are subjected between the time they are laid in the 
fall and hatch in the spring. In the current version, it also implicitly 
accounts for the mortality between hatching and establishment of first 
instars. 



Occasion/ time structure 



Occasion Description 

Initiation Bud burst and shoot elongation complete. 

1 Egg hatch, establishment on foliage. 

2 First ins tars: feeding, growth, and mortality. 

3 Second instars: feeding, growth, and mortality. 

4 Third instars: feeding, growth, and mortality. 

5 Fourth instars: feeding, growth, and mortality. 

6 Fifth instars: feeding, growth, and mortality. Male 
pupation occurs after f ifth-instar feeding. 

7 Sixth instars (female): feeding, growth, mortality, and 
female pupation. 

8 Pupal mortality and adult female emergence. 

9 Eggs laid, aging of foliage, and new foliage potential set. 
10 Redistribution of eggs.!/ Overwinter mortality of foliage 

and overwinter mortality of eggs. 

In Version 3.1, the events of occasion 10 take place in module SQ. 



Table 3 G function assignment by subsystems for the three model versions 



Subsystems 



1. Foliage dynamics 

2. Population dynamics 

3. Growth and feeding 



Subsystems G Functions 



S 2 : 



Versions 1.2 & 2.2 and 
Version 3.1 Lower Module 

1, 4 

2, 5 
3 

6, 7, 8 

9 

10 

11 



12 
13 
14 
15, 16, 17, 18 



19 



Version 3.1 Upper Module 

1 
2 



Identity 



New foliage destruction and 

production. 

Old foliage destruction and (for 

VI and V2) mortality. 

Aging of new foliage into old. 

Determination of daily larval 

mortality. 

Index of days spent feeding on 

old foliage. 

Fecundity or special occasion 

mortality. 

Total mortality (or egg production) . 

Nominal larval growth rate. 
Achieved larval growth. 
Fecundity stress index. 
Food demand and function for 
allocation of demand to new and 
old foliage. 

Index of day on which foliage 
was depleted by feeding. 



Phase or year of the outbreak. 
Redistribution of insects over 
the stand. 

Old foliage mortality. 
Overwinter mortality of eggs. 



Conceptual resolution 



Outbreak, stand 
Phase, tree class 



Phase, tree class 
Occasion, branch 




Time and population count 



TMAX= 1 

where TMAX is the numbei 

outbreaks to be simulated 

; k'{0,... f q-1} 
where q is the number of 
phases, or years, in an 
outbreak. p = <variable>; 
k"6{0,..,p} where p is tt 
number of tree classes in tl 
classification structure. 



Module S| is accessed each 
time k" or k' of module So 
is incremented 



qIO;k fO,... f q-1} 
where q is the number of 
occasions in a phase. 



*ZCOMP sets the values in S t of bj , . . . , b 6 , b 34 , b 35> b 36 , and r t , . . . , r, 2 according to 
the value of z 3 , and rj 3 ,...,r,Q according to the value of 24. Also, the values of 
*5jt X 6jt and X 7j ^ ^o become the initial conditions for x 1t x 2 , and x 3 of Si, respectively, 
and x 4 of S<j is initialized to zero. 



external file 



Figure 3 . --Stand-outbreak model. Version 3.2: Module 
stacking and sequencing. Note: The ovals marked ZCOMP 
and YCOMP represent the FLEX input and output variable 
processor subroutines^ respectively > for each module. 
This figure follows the convention of figures 5 and 6 
in White and Overton (197?). The F's in boxes are 
external data files used in data storage and module 
coupling. For details, see page 46. Couplings between 
SQ and Si are given in terms of variables of SQ. 



FLEXFORMS: VI .2; V2.2; V3.1 MODULE S ; V3.1, MODULE 

FLEXFORM: VI. 2 TITLE: DFTM Population-Branch Model, 

Page 1 of 12 Version 1.2 

March 8, 1974 INVESTIGATORS : Over ton, White, Hunt, 
DATE: 

March 15, 1977 (revised) Colbert 

TEMPORAL RESOLUTION: day, occasion; duration of run, one season. 
SPATIAL RESOLUTION: the hypothetical model branch. 
QUANTITIES MODELED: Biomass (dry weight) of foliage and larva, 

numbers of insects. 

Note: Conceptually this model can be treated as one occasion in an 
annual resolution model, with the annual quantities: 

1. Annual variables 

(i) actual foliage complement 

a) new foliage biomass 

b) old foliage biomass 

(ii) viable egg density 

2. "Phase-specific" parameters 
(i) egg-mass size 

(ii) disease, parasite/predator mortality rates 
(iii) overwinter egg and larval establishment mortality ral 

3. Parameters that originate as "branch-outbreak 11 variables < 
parameters 

(i) nominal foliage complements 

(ii) Phase I initial viable egg density 

(iii) host specific mortality rates 

(iv) (host species) 



FLEXFOKM: VI. 2 
Page 2 of 12 

SUBSYSTEMS: 

1. Foliage dynamics 

2. Population dynamics 

3. Growth and feeding 
DIAGRAM: 



DFTM Population-Branch Model 
March 15, 1977 

Symbol Key: 

O Subsystem designatioi 

O System, subsystem, p] 
function. Numbers ii 
circles are g functic 

I I State variables. 
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VARIABLES AND FUNCTIONS 

1. X List 

Variable Description Units 

X;L New foliage biomass (0 <^ k 62) and g 

new foliage potential for the next 
year (k = 63, 64). 

X2 Old foliage biomass. g 

X3 Viable eggs (k = 0). numbe 

Established first instars (k = 1) . 
Feeding larvae, male and female (1 < k 50). 
Feeding female larvae (51 j< k 60). 
Female pupae (k = 61) . 
Adult females emerged (k = 62). 
Viable eggs produced (k = 63) . 
Viable eggs after overwinter mortality (k = 64). 

x^ Individual larval biomass. mg 

X Individual larval biomass at k = 36. nig 



2. Z Functions 

No inputs other than initial conditions are needed in this 
single-year version of the model. Phase-specific "inputs" are in- 
corporated as parameters. 



3. G Functions 

Functional Form and Description Units 



I" 


[ b 37 . b pi' 


hx 2 +b 38 . ~L . fn 


100 "38 1 




2b QQ 8 19 b 32 f r k C " 
Jo 1 

- 1 g 

otherwise. 



This function accounts for the capacity to produce new foliage 
(g-dryweight) after the completion of feeding by the DFTM (k = 62) and 
before the overwinter needle drop. If feeding results in total de- 
foliation, then new foliage is reduced in proportion to the 
number of days of total defoliation 
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ictional Form and Description Units 





rnn ^ r \ r\ -. r T_ /-i \ 1 -Pj-i-v 1j- A O 


max|0, x 2 - b 38 (l - 1QO ) J for k - 63, 


^ 


g 




otherwise. 



idle drop or litter production (g-dryweight) from the branch. This 
>cess is required to keep the branch in equilibrium under normal 
iditions. 



x- for k = 62, Aging of new foliage to old 

foliage. 



otherwise. 



min(g_,, x n ) for 1 < k < 60, 
ID 1 

g 
otherwise. 

struction of new foliage (g-dryweight) for each day of feeding (1 <_ k 

min(g 1Q - g , x ) for 1 < k < 60 and g < g , 
lo 4 Z 4 ID 

g 
otherwise. 

truction of old foliage (g-dryweight) for each day of feeding (1 k 
which demand for food can not be fully satisfied by feeding 
new foliage. 

(1 - r ) (1 - r ) (1 - b c ) for 1 k 1 60, 

S l S l b S l .-1 

d 

otherwise. 

survivorship per day (d ) associated with disease mortality (r ) , 

S l 
dator and parasite mortality (r ,), and background mortality 

S i "T" D 

). Each is instar specific (see special function S..). 
1 
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1 - r 







12 



for 1 1 k 60, 
otherwise. 



Dhe daily survivorship associated with stress from feeding 
m old foliage; instar specific. 



j = f 8 6 8 7 15 ~ X 

u 

Negative total daily mortality rate. 


for 1 k 60, 
otherwise. 
[Assume = ! ] 


1 -jc* + P- 
l x l ^ 817 


for k = or k - 63, 
otherwise. 



(implicitly, x* (k+1) 



x* (k) + g 1? (k)) 



Cotal number of days (d) spent feeding on old foliage during 
demands are fully satisfied. 



o 



g 8 ) 



g 14 (i - 



-b 



35 







for k = 0, 
for k = 50, 
for k = 61, 
for k = 62, 
for k = 63, 
otherwise. 



This function accounts for changes in number of larvae from 
establishment mortality (k = 0) , male pupation (k = 50), female 
pupal mortality (k = 61) , viable eggs produced by the remaining 
females (k = 62) , and overwinter mortality (k = 63) . Fecundity is 
reduced in proportion to the number of days spent feeding on old 
foliage (g~) and the final weight of the female larva (g-,/)- 
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Functional Form and Description 



=11 



Larval mortality, male pupation, 
pupal mortality, egg production. 



'12 



for 1 k <_ 35, 
for 36 <_ k 60, 



for x < or 
x + -x. <_ or 
k = or k > 61. 



The nominal instantaneous 
larval growth rate. Lar- 
vae grow at the nominal 
rate as long as food is 
available. 



13 



Here 



25 



for k = 0, 

for 1 k ^ 60, 

for k > 61. 



for (g 



g !5 g !6 + g 17 g !8 



15 



0, 



g 



18 



for (g 15 + g 1? ) = 0, 



Achieved larval growth increment. The proportion <J> is the 
ratio of actual foliage destruction to total foliage demand; 
if less than one, then the larvae did not feed for the full 
day and hence the growth is reduced according to the propor- 
tion <j>. 



mg 



max 



24 



Fecundity stress index, none 
which is measured relative 

to the nominal final biomass(.b^, ) . 

24 



1 if 8 16 < x x 

and 1 <_ k <_ 60, 
otherwise. 



The indicator of a day spent 
feeding exclusively on new d 
foliage. 
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fieri 

3 1 +6 L b 21 



f , 

(b 



20, , g !2 

T"~> (e 

g !2 



16 = 1 for 1 k <_ 60, 

otherwise. 

New foliage demand for the day's feeding and growth. 



'17 



i if g 15 = o, 

1 k 60, and (x + 
otherwise. 



> g 



lg , 



The indicator of a day spent feeding on old foliage with 
demand completely satisfied. 



10 



,. 
(b 1 



8 



12 



S;L + iz I & 21 - 19 S 12 3 4 

for 1 _< k _< 60, 

otherwise. 

Old foliage demand for the day's feeding and growth. 



60 - [x* + g 9 ] for k = 62, 



otherwise. 



(implicitly, x* (k+1) = x^ (k) -f g 15 (k)) 

Total number of days during the larval period for which food 

demands were not fully satisfied. 
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for k = 0, x 3 0, x + x 2 0, or k _> 61, [1] 



4/ 



for k = 1 or [x 4 >. (1 - 



n(x 4 /b 25 ) 
k 1 ! 



^g^ and 2 k 35], [2] 
- b 31 )g 21 and 2 k 35, [3] 



b 23 for k = 36 or [x 4 >_ (1 - b 31 )g 21 and 37 k 60], [4] 



n(x 4 /x 5 ) 
L. k-36 



)g and 37 k 60. [5] 



; modified growth rate of an individual larva. [3] and 
| are growth rates for trajectories below the normal; [2] 
1 [4] for normal trajectories. 



for < k < 35, 



35b + (k - 35)b 
b_,e for 36 < k < 60, 





otherwise. 



ninal larval growth trajectory. 



4/ 



12 



i: 



if k = 36, 



otherwise. 



Larval biomass after 35 
days of feeding, bench- 
mark for determining re- 
duced growth rate. 



4/ 



= max 



f P4 + g 13~| 

2 g " lf 

[ L S 21 J 



Fecundity stress index. 



- Note: These four g-functions are not operational, but are given here 
exhibit how variable growth rates for the larvae might be implemented if 
itional food becomes available after foliage is depleted. The functions g 
would replace the present g,- an< 3- 814* respectively. 
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Units 





f l,l g l g 4^ 


g 




f l 2 = 3 


g 




f 2 2 = -e - g 

^ 5 * St: oo 


See state variables X list, for g 




f 3 3 = g ll 

O,O J_J- 


descriptions. number 




f = e 
4,4 g !3 


mg 




f 5,5 * g 22 j 


mg 




5. 


Y Functions 


Functional Form and Description Units 



y . = x. , 
'i i 



i = 1, 2, 3, 4 



See X list for descriptions 
and units. 



6. Special Functions 
Functional Form and Description 



Units 



s r = 



10 



The current ins tar of the 
population. 



none 



PARAMETERS 



7. B Parameters 



identical to 






of Version 2.2. See 



page 7, 8 of Version 2.2 FLEXFORM. 



Note: For this function, [x] is the largest integer that does not 
exceed the magnitude of x. 
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List 



23 



24 



'25 






32 



33 



34 



Value 
0.0625 

52.063 
0.197 



b I, .50; II, 0.62; 
III, .75; IV, .90 

b I, 200; II, 200; 

III, 150; IV, 150. 

bo (not used) 
b (not used) 

(not used) 



Description 

Instantaneous growth rate for 
36 <_ k 60. 

Nominal biomass of female larvae 
at the end of the sixth instar, 
k = 61. 

Initial nominal biomass of larvae 
after establishment. 

Phase-specific pupal mortality 
rate. 

Phase-specific nominal egg- 
mass size. 



1.0 X 10 



0.02 



0.01 



0.50 



b I, .50; II, .60; 

III, .85; IV, .90. 



1, 2, 3 or 4 



6/ 
""" 

6/ 
~ 



38 



39 



0.0 



Error tolerance in larval 
growth traj ectory . 

Coefficient of reduction in 
bud capacity from total branch 
defoliation. 

Coefficient of reduction in 
fecundity from feeding on old 
foliage. 

Proportion of larvae that are male 
at the end of the fifth instar. 

Overwinter mortality. 

Phase index. 

Nominal percent of new foliage. 

Nominal total foliage biomass 
per model branch. 

Phase-specific larval establish- 
ment mortality. 



mg 



The values of b^j and b^g are chosen according to the properties c 
host tree for a particular simulation. 
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8. 


R Parameters 


List 


Value 


I 


II 


III 


IV 


Ins tar 


r 


1 


0. 








.0 


0. 


002 


0. 


025 


Ist^ 


r 


2 


0. 








.0 


0. 


003 


0. 


028 


2nd 


r 


3 


0. 








.0 


0. 


006 


0. 


031 


3rd 


r 


4 


0. 








.0 


0. 


013 


0. 


034 


4th 


r 


5 


0. 








.001 


0. 


035 


0. 


035 


5th 


r 


6 


0. 








.001 


0. 


028 


0. 


035 


6th 

* 


r 


7 


0. 








.0 


0. 


001 


0. 


005 


1st' 


r 


8 


0.0 





.0 


0. 


002 


0. 


006 


2nd 


r 


9 


0. 





0.0 


0. 


003 


0. 


007 


3rd 


r 




0. 





0.0 


0. 


010 


0. 


021 


4th 




10 




















r 




0. 








.001 


0.016 


0. 


033 


5th 




11 




















r 




0. 








.001 


0. 


042 


0. 


056 


6th 




12 


















j 


Douglas-fir 


Grand 


fir 


r 






1. 









1.0 






1st" 




13 




















r 






0. 


75 






0.8 






2nd 




14 




















r 






0. 


10 






0.5 






3rd 




15 




















r 






0. 









0.1 






4th 




16 




















r 






0. 









0.0 






5th 




17 




















r 


1 O 




0. 









0.0 






6th 



Description 



Ins tar and phase- 
specific daily 
disease mortality 
rate. 



Ins tar- and phase- 
specific daily 
predator /parasite 
mortality rate. 



Instar- and host- 
specific daily 
food stress 
mortality rate. 
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II 



Initial Conditions' 



x l 


50.4 


X 2 


93.6 


X 3 


3.5 


X 4 


0.0 





CNERAL RUN INFORMATION 

TSTART = 0_ 

TMAX = 64- 

q =1 

LP, DUMP-( All y's, all g's 

LP, DUMP frequency: JL 

RUN LOG: This version of the model was developed for documentati 
purposes and was run only to verify model output 
against V3.1. 

Source code file: DAILY 

9/ 
FLEX3 overlay- : *DAILY 



The full range of admissible initial conditions is available from th 
thors upon request. 

-/ LP stands for "line printer". DUMP stands for "dump file", a file us 
r compact storage of the details of a simulation for later review or recon 
ruction. 

9/ 

"Overlay" as used here is the machine-executable code, and is called 

bsolute element" on some other computers. 



FLEXFORM: V2.2 TITLE: DFTM Population-Branch Model 

Page 1 of 11 Version 2.2 

INVESTIGATORS : Over ton, Colbert, Whit 
DATE: December 1, 1975 

March 15, 1977 (revised) 

TEMPORAL RESOLUTION: instar, occasion; duration of run is one seasc 
SPATIAL RESOLUTION: the hypothetical model branch 
QUANTITIES MODELED: biomass (dry weight) of foliage and larvae; 

number of insects 
SUB S Y STEMS : 1 . Foliage dynamics 

2. Population dynamics 

3. Growth and feeding 

DIAGRAM: 

(See the diagram on page 2 of Version 1.2 FLEXFORM.) 

VARIABLES AND FUNCTIONS 

1. X List 

List Description Units 

x- New foliage biomass g 

x ? Old foliage biomass g 

x Viable eggs, k = 0. number 

Established first instars (k = 1). 

Feeding larvae, male and female (k = 2, 3, 4, 5). 

Feeding female larvae (k = 6). 

Female pupae (k = 7). 

Adult female moths (k = 8). 

Viable eggs produced (k = 9). 

Viable eggs after overwinter mortality (k = 10). 
x. Individual larval biomass mg 
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2. Z Functions 

No inputs, other than initial conditions, are in this sir 
year version of the model. Phase- specific "inputs" are incorp 
as parameters. 



3. G Functions 

Functional Form and Description- 

V^ Pi + *a + " . . 1 for fc - 8, 



100 



Pi + 
L 



otherwise. 

New foliage capacity. 



(max 1 , 




x - b ^(1 - in ) I for k = 9, Needle drop. 

otherwise. 



- P 

" lo 



x for k = 8, Aging of new foliage. 



g 3 

1 ^ otherwise. 



,, x ) for 1 <_ k _< 6, Destruction of 
\ new foliage, 

otherwise. 



c 



, x ) for 1 k 6, Destruction of ol 

foliage. 

otherwise. 



' For more complete descriptions see Version 1.2 FLEXFORM. 
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Units 



S 6 = 



(1 - r k ) (1 - r k+6 ) (1 - b fc ) for 1 k 6, 







otherwise. 



Daily survivorship after disease, predator /parasite, and background 
mortality. 



* P +12 



otherwise. 
Daily survivorship after food stress. 



10 10- g 



15 



- 1 for 1 k 6, 
otherwise. 







Negative total mortality rate for an instar. [Assume = 1. ] 



r 

lx* + 



for k = or 9, 
otherwise. 



(Implicitly, x*(k + 1) = x*(k) -1- g 1? (k)) 

Total number of days spent feeding on old foliage during 
which demand is fully satisfied. 



36 



Special Occasion Processes: 
for k = 0, Establishment mortality. 



-b n (l + g Q ) for k = 5, Male pupation. 

JU o 



-b 



34 



for k = 7, Female pupal mortality. 



b 35 8 14 (1 - b 28 8 9 ) ~ 



-b 



37 



for k = 8, Viable eggs produced, 
for k = 9, Overwinter mortality, 
otherwise. 
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g = (g ft 4- g )x Larval mortality, male pupation, 
pupal mortality, egg production, 
and overwinter mortality. 


numbe 



'12 



for 



<_ 0, 



~ 



0, ork>6, Instantaneous 



22 



for 1 <_ k < 4, 
for k = 4, 



daily growth 
rate. 



b 0/ for 4 < k < 6. 
^ 24 





r b 25 for 


k = 0, 




g 13 = 


[10<j)g ] 
x 4 (e - 1) for 


1 k 5 6, 


mg 




for 


k > 6, 






where 



Achieved larval growth increment for an ins tar. The pro- 
portion cf> is the ratio of actual needle destruction to 
total needle demand; if less than one, the larvae did not 
feed all 10 days, and the number of days spent growing is 
reduced according to the proportion cj) . 



'14 



max 




Fecundity stress index. 



'15 



where a e (0, 1, ...,10} is 

the largest integer such that 

UVl(k, x, b, a) < x , 

for k =~0 or k > 6 or if no a exists that 

satisfies the above inequity. 



Time spent feeding exclusively on new foliage. 
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Unit; 



UVl(k, x, b, min[g 15 + 1, 10]) 



for x > and 1 k 6, 
otherwise , 



Total new foliage demand. 



10 - g 




15 



where n e {0, 1, . . . , 10 - -,,-} 

is the smallest integer such that 

[UV2(k, x, b, g 15 , n + 1) - 

( X;L - UVl(k, x, b, g 15 ))l > * 2 and 

otherwise for 1 <^ k _<_ 6, 

for k = or k > 6. 



The number of days spent feeding on old foliage for which demand is 
fully satisfied. 



g 



18 



max[UV2(k, x, b, g^, 10 - g^) - 

(^ - UVl(k, x, b, 8 15 )), for 1 _< k 6, 



for k = or k > 6. 

Total old foliage demand for remainder of the ins tar. 

for k = or 9, 
for k = 7 or 8, 
10 - [g, r + g, -,] otherwise. 



X* 

X 2 



(Implicitly, x*(k + 1) = x*(k) + 10 - [g 15 (k) + g iy (k)]). 

Total number of days during which demand for food is not 
fully satisfied. 
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4. 



F Functions 



Functional Form and Description 



Un 



f l,l s l g 4 




f l,2 = g 3 


See state variables, X list, 


f = _0 _ 

2,2 S 5 S 2 


> for descriptions. 


= 2 


num 


3,3 11 




f 4,4 - g !3 


m 





5. Y Functions 



i - v 



i = 1, 2, 3, 4 See X list for description 
and units. 



6. Special Functions 
Functional Form and Description 



21 



= b 



6+k 



r g i^i a 

- Lv J , 



i ~ 



and 



Un 



UVl(k, x, b, a) = U(x, b, k) . ^(b, a, k) , 

where o T. , 

10" 3 f b 2ol g !2 

' b ' k) - Sr b i9 + z: < e -"^ 



UV1 = for k i { 1, ..., 6}. 

Total new foliage demand for an ins tar. U is the food demand 
for metabolism and growth, V]_ is the product of the destruction/ 
consumption ratio and the convolution of daily larval growth and 
daily larval survivorship during feeding on new foliage. See 

14 and 
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Functional Form and Description 



V 2 (b, a, n, k) = 



a ag 



12 



1 - 


g !2 
g 6 g ? e _ 


n 


1 - 


g !2 
f6 g 7 e . 





Units 



UV2(k, x, b, a, n) = U(x, b, k) . V (b, a, n, k) . 

~ ~ ~ z. ~ 

where U is given in UV1 above 
and 



UV2 = for k i { 1, ..., 6} 



Total old foliage demand for an ins tar. U is the food demand 
for metabolism and growth, ^2 ^ s tne product of the destruction/ 
consumption ratio with the convolution of the daily larval growth 
and daily larval survivorship during feeding on old foliage. 



PARAMETERS 



7. B Parameters 
List Value 



!2 



14 
15 
16 

17 
!8 



0.02 
0.02 
0.02 
0.02 
0.02 
0.02 

5.4 

6.25 

6.25 

2.71 

2.27 

2.0 

5.4 

6.25 

6.25 



69 
29 



3.0 



Instar 



1st 
2nd 
3rd 
4th 
5th 
6th 

1st' 

2nd 

3rd 

4th 

5th 

6th. 

1st' 

2nd 

3rd 

4th 

5th 

6th. 



Description 



Unit 



Instar-specif ic background 
daily mortality rate. 



Instar-specif ic 
destruction /consumption 
ratio for feeding on new 
foliage. 



Ins tar-specif ic 

destruction/ consumption ratio 

for feeding on old foliage. 



non< 



noi 



j. age. u u 


J- J. JL 


rid. i. uu J.J 5 JLZ7 I l 




List 


Value 


Description 


Unii 


b !9 


1.19 


Assimilated food, growth coefficient. 


noi 


b 20 


0.081 


Assimilated food, respiration coef- 


cf ' 


A.VJ 




ficient. 




b 21 


0.10 


Assimilation efficiency. 


noi 


b 22 


0.1147 


Instantaneous growth rate for first 


d~ " 


<. <. 




three ins tars. 




b 23 


0.0886 


Instantaneous growth rate for fourth 


d" 1 






instar. 




b 24 


0.0625 


Instantaneous growth rate for fifth and 


d" 1 






sixth instars. 




b 25 


0.197 


Initial biomass of larvae at hatch. 


mg 


b 26 


52.063 


Nominal biomass of female larvae 


mg 


z.u 




at the end of the sixth instar. 




b~_ 


0.02 


Coefficient of reduction in bud 


d 


27 




capacity from total branch defoliation. 




b 28 


0.01 


Coefficient of reduction in fecundity 


d 






from feeding on old foliage. 




b 29 


1,2,3,4 


Phase index. 


not 


b 30 


0.50 


Proportion of larvae that are male 


not 


-J W 




at the end of the fifth instar. 




H/ 












Nominal percent new foliage. 


noi 


b 31 ll/ 












Nominal total foliage biomass. 


g 


b 33 


(not used) 






b 34 


1:0.50; 


Phase-specific pupal mortality rate. 


occ< 




11:0.62; 








111:0.75; 








IV : . 80 






b 35 


1,11:200; 


Phase-specific maximum egg-mass size. 


nur 




III, IV: 150. 






b 36 


I, II, III, 


Phase-specific establishment mortality 


OCCc 




IV:0.0 


rate. 




b 37 


1:0.50; 


Phase-specific overwinter mortality 


OCCc 




11:0.60; 


rate. 






111:0.85; 








IV:0.90 









Note: bji and b^2 have values that depend on the host tree being 
simulated and are chosen, accordingly, prior to a specific simulation. 
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List 


Value 


r i 


I 
.0 


- Phase - 
II III 

0.0 0.002 


IV 
0.025 


Insta: 
1st" 


I 




















r ? o 


.0 





.0 





.003 




0. 


028 


2nd 


r!r 


.0 





.0 





.006 




0. 


031 


3rd 


t 


.0 





.0 





.013 




0. 


034 


4th 


r s 


.0 





.001 





.035 




0. 


035 


5th 


C 6 


.0 





.001 





.028 




0. 


035 


6th. 


c 


.0 





.0 





.001 




0. 


005 


1st" 


c 


.0 





.0 





.002 




0. 


006 


2nd 


C Q 


.0 





.0 





.003 




0. 


007 


3rd 


r 


.0 





.0 





.010 




0. 


021 


4th 


: v? 


.0 





.001 





.016 




0. 


033 


5th 


1 1 




















: 12 


.0 





.001 





.042 




0. 


056 


6th^ 


Douglas-fir 


grand 


fir 








1 


.0 






1 


.0 




1st" 


: 







.75 









.8 




2nd 


.14 







.10 









.5 




3rd 


15 



























.0 









.1 




4th 


: 17 







.0 









.0 




5th 


: ' 







.0 









.0 




6th , 



Description 



Units 



Instar- and phase- 
specific daily 
disease mortality rate 



Instar- and phase- 
specific daily 
predator /parasite 
mortality rate 



,-1 



Instar- and host- 
specific daily food- 
stress mortality rate. 
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9. Initial Conditions 


X- 


50.4 


1 




x n 


93.6 


2 




X "5 


75.0 


3 




x, 


0.0 


4 







GENERAL RUN INFORMATION 

TSTART = 0^ 

TMAX = 10 

q = I 

LP, DUMP: All y's; all g's 

LP, DUMP frequency: 1. 

RUN LOG: This version of the model was developed for the 

purpose of documentation and was run only to verify 
this model form against V3.1. See appendix B for v< 
fication output. 

Source code file: INSTAR 
FLEX3 overlay: *INSTAR 



These initial conditions, along with the phase II b and r param 
for grand fir, were used to develop the appendix B verification output. 



FLEXFORM: V3.1 TITLE: DFTM Stand- Outbreak Model Versior 

MODULE : S 

Page 1 of 7 Outbreak/ Stand Module. 

INVESTIGATORS : Qverton, Colbert, White 
DATE: September 17, 1976 

March 15, 1977 (revised) 

TEMPORAL RESOLUTION: Outbreak/phase. 

SPATIAL RESOLUTION: Stand (by tree class). 

QUANTITIES MODELED: Biomass (dry weight) of foliage and larvae, numbei 

of insects per model branch, for each tree class. 

DIAGRAM: See figure 3. 

VARIABLES AND FUNCTIONS 
1. X List 

List Description Units 

Phase 



x 2 II 

x 3 III 

x. IV 
4 



Weighted mean number of eggs number 
per model branch over the tree 
classes of the stand. 



The vector variables below contain an element for each tree class: 

x Actual new foliage potential (on entry to S.. ) . g 

Actual new foliage biomass (on exit from S ). 

x. Actual total foliage biomass (on entry to S ) . g 

Actual old foliage biomass (on exit from S ). 

x 7 Number of viable eggs. - number 
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DULE: S Q 
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iDULE: S March 15, 1977 



13/ 
Z Functions 

st 

The number of tree classes in the stand for the 
current model run. This number is entered by 
the user during a model run. This value is stored 
in b57 for use in computing the final weighted mean 
number of eggs. This is also the value of p in the 
input command file and must be set before a model run. 

The vector variables below contain an element for each 
tree class: 

Nominal percent new foliage. 

3 Nominal total foliage biomass. 

Number of trees. 
4 

Weighting factor. 

f Host species (1 = Douglas-fir, 2 = grand fir). 
6 



G Functions 



nctional Form and Description 



= k + 1 Phase index. 



14/ , f X 7j 

3 ~L 



b 6 2 (x k f ~ X 7j) for k? 



otherwise. 

distribution of viable eggs over the stand. (Implemented by 
call to special function REDIST.) 



'^ Note: the z values are entered by the user in specifying a stand- 
:break model run. 

147 

Functions 2, g3 , and 4 process all tree classes on the first pass 

rough the g functions each phase. 



3DULE: S Q March 15, 1977 

ige 3 of 7 

mctional Form and Description Unit; 

for k' > 0, 



otherwise. 

Ld foliage survivorship. 

)g 2 . for k' > 0, 

numb< 
[0 otherwise. 

/erwinter survival of eggs. 

Y Functions 

YCOMP calls REDIST for computation of the mean fall egg density 
/er the stand after phase IV. No y functions, per se, are included in 
lis module. Instead, YCOMP accesses certain special functions to produ 
)ecially formatted output tables. These are documented elsewhere 
see: Logical Unit Manipulation by User-Specified Functions). 

Special Functions 
ibroutine REDIST (k ! , n, b, x) 

This subroutine first computes the mean egg density over 
le stand as follows: 

v = i 

lere 

N 

T = y z/ .z_. 

L- 4j 5j 

id N is the number of tree classes; z^. is the number of trees in the 
:h tree class; z,. . is the weight associated with a tree of the j th tree 
.ass; x is the number of viable eggs laid for the jth tree class; and 
., is the fall mean egg density over the stand in phase k* . 



The subroutine then redistributes the eggs over the stand according 



2' 



MODULE: S Q 
Page 4 of 7 

PARAMETERS 


6. B 


Parameters 




Value 


b, 


0.0 


i 




t_ -'- 


0.0 


9 




t>0 


0.0 


X 




b 4 


0.0 


b- 


0.0 


s 




b* 


0.0 


6 




b- 


0.0 


7 




b fi 


0.0 


"h 


0.0 


b 


0.0 


b* 1 

b !2 


0.001 
0.001 


b !3 


0.002 


J 14 
t 15 
h 16 

D-. _. 


0.003 
0.006 
0.013 
0.035 


i, 17 
b ]8 


0.028 


b 19 
b 20 

hi 

b 22 

b 23 
b 24 


0.025 
0.028 
0.031 
0.034 
0.035 
0.035 


b 25 
h 26 
b 27 
b 28 

b 29 
b 30 


0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


b 31 
b 32 
b 33 

^ 


0.0 
0.0 
0.0 
0.0 
0.001 


b 


0.001 



March 15, 1977 



Description 



Ins tar 

Istl 

2nd 
3rd 
4th 
5th 
6thJ 



Phase 



2nd 
3rd 
4th 
5th 
6th ^ 

IstT 

2nd 

3rd 

4th 

5th 



II 



III 



1st" 

2nd 

3rd 

4th 

5th 



IV 



Ins tar-specif ic 
disease mortality 
rate. 

(r 1? . . . , rg in 
Version 2.2) 



2nd 
3rd 
4th 
5th 



1st" 

2nd 

3rd 

4th 

5th 

6th 



II 



Ins tar-specif ic 
predator /parasite 
mortality rates. 
(r 7 , . . . , r 12 in 
Version 2.2) 



(continued) 
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Value 



DFTM Stand-Outbreak Model 
March 15, 1977 



Description 



Units 



Ins tar 



Phase 



0.001 


1st" 


^ 


0.002 


2nd 




0.003 


3rd 




0.010 


4th 


in 


0.016 


5th 




0.042 


6th^ 




0.005 


1st' 




0.006 


2nd 




0.007 


3rd 




0.021 


4th 


IV 


0.033 


5th 




0.056 


6th^ 


J 


0.50 


I 


0.62 


II 


0.75 


III ' 


0.80 


IV . 


200 


f 


200 


II 


150 


III ' 


150 


IV 



Pupal mortality rate, 
(b , in Version 2.2) 



Maximum egg-mass size 
eggs /mass 0^35 in 
Version 2.2) 



Number of tree classes (see z functions) 



occasion 



number 



number 



Phase 



0.50 
0.60 
0.85 
0.90 

0.0 

0.0 
0.0 
0.0 
0.0 



I 

II 

III 

IV 



Overwinter 
mortality rate. 
(b 37 in Version 2.2) 

Redistribution dispersal factor 



I 

II 

III 

IV 



Establishment 

mortality rate. 

(b 36 in Version 2.2) 



occasion 



none 



occasion 
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7. R Parameters 





Value 


Description 


Un 






Instar 




r 


1.0 


1st 1 




r 9 


0.75 


2nd 


Instar-specif ic daily food 


n 


L 


0.10 


3rd 


stressi^./ mortality rates for 




* L 


0.0 


4th 


Douglas-fir host species. 




H- 

r 5 


0.0 


5th 


^ r !3 > r !8 *- n Version 2.2) 




r 6 


0.0 


6th^ 






r_ 


1.0 


Ist^ 






7 


0.8 


2nd 


Instar-specif ic daily food 


n 


r q 


0.5 


3rd 


stressi^./ mortality rates for 




r io 


0.1 


4th 


the grand fir host species. 




JLU 
r ll 


0.0 
0.0 


5th 


(r^3, ..., r-^g in Version 2.2) 




r !2 




* 




r 


0.02 


1st" 






13 

r l4 


0.02 


2nd 


Instar-specif ic daily background 


n 


JLH- 

r !5 


0.02 


3rd 


mortality rate. (b^ > bg in 




JL _J 

r !6 


0.02 


4th 


Version 2.2) 




JLVJ 


0.02 


5th 






r !7 
r !8 


0.02 










8. Initial Conditions 



See appendix C for an example of the specification of initial 
conditions. 



' Stress from feeding on old foliage during early instars. 
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GENERAL RUN INFORMATION 

Module Sequence Number = 

TMAX = !_ q = 4. 

Module type: ghost 

Variables to be monitored: none. 

Variables to be dump filed: none. 

FILE NAMES 

FLEX2 function overlay name: *V31SO 
Source file name: V31SO 

COMMENTS 

Model behavior will be documented in a separate paper. A record 
the model runs made in studying behavior will be included therein. 
Appendix A contains the source code for Version 3.1 and the command f 
for initializing a simulation. Appendix C contains the inputs and OL 
from a verification simulation. 



FLEXFORM: V3.1 TITLE: DFTM Stand- Outbreak Model, 
MODULE : S Version 3.1, Tree/Branch Module 

Page 1 of 5 INVESTIGATORS: Overtoil, Colbert, White 

DATE : September 17 , 1976 

March 15, 1977 (revised) 

TEMPORAL RESOLUTION: Phase/instar , occasion. 
SPATIAL RESOLUTION: Tree class/hypothetical model branch. 
QUANTITIES MODELED: Biomass (dry weight) of foliage and larvae, 

of insects. 
SUBSYSTEMS: 1. Foliage dynamics 

2. Population dynamics 

3. Growth and feeding 
DIAGRAM: 

(See the diagram on page 2 of Version 1.2 FLEXFORM.) 

VARIABLES AND FUNCTIONS 

1. X List 

List Description Units 



x New foliage biomass (k 1 =0, ..., 8). 

New foliage potential for next year (k 1 = 9, 10). 

x Old foliage biomass. 

x Viable eggs (k' = 0). 

Larvae (k' = 1, . . . , 5) . 
Female larvae (k f = 6). 
Female pupae (k' = 7). 
Adult female moths (k 1 = 8). 
Viable eggs produced (k 1 = 9, 10). 

x. Individual larval biomass. mg 
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2. Input Variables 
Z Functions 
Functional Form 



March 15, 1977 



Z l = Z 2j 



in S 




Q 



in S 



Z 4 = Z 6j 







Description 

Nominal percent new foliage. 
Nominal total foliage biomass, 
Phase index. 
Host species index. 



Units 

none 

g 
none 

none 



Parameter Assignments 
B Parameters 



respectively, for z = 1, 2, 3, 4. 



b 49 in S f 2 3 = 



'34 



50 ln S f Z 3 = 2 ' 

51 in S Q if z 3 = 3, 

52 in S Q if z 3 - 4. 



Value, Description, and 

See V3.1 FLEXFORM, Modul 
S n , for the values, desc 
tions, and units of the 
assigned b and r paramet 



35 



b 53 in S Q if z 3 = 1, 
b 54 in S " Z 3 = 2 ' 



in S Q if 



3, 



b 56 in S " Z 3 = 4 ' 



36 



b 63 n S if 2 3 = ls 



in S Q if 



2, 



b 65 in S Q if z 3 = 3, 
b 66 ln S if Z 3 = 4 " 
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R. Parameters 





'[b 


1' ' 


.., b 6 


i 


in 


S 


0'. 


respectively, 


for 


Z 3 


= 1 


[r , ..., r ] - < 


[b 


7' ' 


' b l 


9 


in 


S 


0' 


respectively, 


for 


Z 3 


= 2 


X o 


[b 


13' 


. . . , b 


18 ] 


in 


S 


0' 


respectively , 


for 


Z 3 


= 3 




[b 


19' 


, b 


24^ 


in 


S 


0' 


respectively, 


for 


Z 3 


= 4 




'[b 


25' 


, b 


30 ] 


in 


S 


0' 


respectively, 


for 


Z 3 


= 1 


[i 


[b 


31' 


. . . , b 


36 - 1 


in 


S 


0' 


respectively, 


for 


Z 3 


= 2 


T" T" = ^ 
J._, ., .L _ ~ J 


[b 


37' 


. . , b 


42 ] 


in 


S 


0' 


respectively, 


for 


Z 3 


"" 3 




Jb 


43' 


. . , b 


48 ] 


in 


S 


0' 


respectively, 


for 


Z 3 


= 4 




[r 


1' 


..., r 6 


] 


in 


V 


respectively, 


for 


Z 4 


= 1 


[r 13 , ..., r 18 ] = 


[r 


7' 


'' r !2 ] 


in 


V 


respectively, 


for 


Z 4 


= 2 





3. 



G Functions 



g Q are the same as in model Version 2.2, except that g^ and 
the portion of g]_Q accounting for occasion 10(k=9) have 
moved to the upper module and the following substitution 
have been made : 



Version 2.2 
b 31 



Version 3.1, S 



32 



"29 "3 

Because processing in this module takes place using the 
lower level, k is replaced by k 1 in all of the g functio 
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1. F Functions / 



Functional Form and Description Units 

- 84) S 

See state variables, g 



-10g q X List, for description g 

10g - by corresponding indexes. number 

10g 13 nig 



. Y Functions 

. = x., i = 1, 2, 3 See X list for description and units. 



D. Special Functions 

The special functions UV1 and UV2 are the same as those used in 
Version 2.2 except that k in the call string has been replaced by k 1 . 



AMMETERS 



For values and definitions, see Version 2.2 FLEXFORM. 
B Parameters 

Description 



- , ..., b. Carried in S n and assigned as a variable input. 

*y 3 . . . , b-i o 

13 ' "'18 As in Version 2.2 FLEXFORM. 



FLEXFORM2 V3.1 
MODULE.* S 
Page 5 of 5 



DFTM Stand-Outbreak Model 
March 15, 1977 



R Parameters 



r ? , ..., 
r 13 , ... 



Carried in S and assigned as a variable input. 



9. Initial Conditions 

See appendix C for an example of the specification of initial 
conditions. 



GENERAL RUN INFORMATION 

Module Sequence Number = 1. 
TMAX =1 q = 10 TARGET 
Module type: terminal. 
Variables to be monitored: y- , 
Variables to be dump filed: y- : 



V y 3- 



Resolution: 
Resolution: 



FILE NAMES 



FLEX2 function overlay name: *V31S1. 
Source File name: V31S1. 



COMMENTS 

Model behavior will be documented in a separate paper. A reco 
of the model runs made in studying behavior will be included there! 
Appendix A contains the source code for version 3.1 and the command 
file for beginning a simulation. Appendix C contains the inputs 
and outputs from a verification simulation. 



LOGICAL UNIT MANIPULATION 
BY USER-SPECIFIED FUNCTIONS 

The DFTM stand-outbreak model, V3.1, makes extensive utilization oJ 
ogical units (LUNs) in storing various vector values and in producing 
pecially formatted output tables. Much of the information concerning 
hese LUNs is summarized in table 4 and figure 4. This section details 
he information stored on some of the LUNs and the user specified 
unctions that produce that information. 



Table 4. Logical units: assignment list 

Note: All LUNs are equipped as linear access files except 
LUN 25, which is a random access file. 



Logical Unit 
Number (LUN) 


Description 


Data Format 


1 


Phase I initial conditions 


(12, 7F11.5, 13) 


2 


Phase I outputs 


(3F12.6) 


3 


Phase I redistributed insects 


(1F12.6) 


4 


Phase II initial conditions 


(12, 7F11.5, 13) 


5 


Phase II outputs 


(3F12.6) 


6 


Phase II redistributed insects 


(IF12.6) 


7 


Phase III initial conditions 


(12, 7F11.5, 13) 


8 


Phase III outputs 


(3F12.6) 


9 


Phase III redistributed insects 


(1F12.6) 


10 


Phase IV initial conditions 


(12, 7F11.5, 13) 


11 


Phase IV outputs 


(3F12.6) 


12 


Phase IV redistributed insects 


(1F12.6) 


13 


Mortality and top-kill for 


(1F7.4, 6F6.3) 




Douglas-fir host species 




14 


Mortality and top-kill for 


(1F7.4, 6F6.3) 




grand fir 




15 


Table 1: Insect-foliage 


variable BCD 




outbreak summary 




16 


Table 2: Foliage-stand 


variable BCD 




outbreak summary 




17 


Table 3: Simulation run model 


variable BCD 




parameterization 




18 


Screen display format of table 1 


variable BCD 


19 


Scratch pad for temporary storage 


(1F12.6) 




of module SQ function g^ output 






values (3.) 




20 


Table 4: ""(optional) detail 


variable BCD 




output of state variables during 






each phase 




21 


State variable values at each 


(4(2X, F14.7)) 




occasion for each tree class and 






phase in a simulation 




22 


Screen display format of 


variable BCD 


23 


table 2: part I 
Screen display format of 


variable BCD 




table 2: part II 




24 


Screen display format of 


variable BCD 




table 2: part II 




25 


State-variable values : input 


binary 



ILUN I? 


3 
s 




SLUM 16 


2 

< 


f LUN 24 


LUN 15 
Lineprinter 
Table 1 


1 LUN 23 


2 
Li 


ILUN 22 


2, 


LUN 18 
Interactive 
terminal 
Table 1 


S> 

x^* 





I LUN 10 


JW. 


I LUN 7 


III 


| LUN 4 


LUN 1 Annual 
resolution : 
Phase specific 
inputs 
Phase I 


II 





^ 


S* 


LUN 25 (RAF 
Same as LUN 


) 
21 
but in 
Binary 


LUN 21 
Detailed state 
variable values 
output in BCD 


LUN 20 
Table 4 

Detail run 
information 


< V 


EXTERNAL\ 
PROGRAM \ 
TABLE4 Y*"""*"" 




\ 


^ 



Figure 4 . --Stand- outbreak model,, Version 3.1, 
LUN assignment and identity. 



LUNs 1, 4, 7, and 10 contain phase-specific initial condition 
information for phases I, II, III, and IV, respectively. The variabl 
identity, by record on the LUN, is: 

Z 6j> Z 4j> Z 5j> Z 2j> Z 3j> X 5j> X 6j> X 7j ' j 

rfiere j is the tree-class identification index, and all variables are 
those of module SQ. 

LUNs 2, 5, 8, and 11 contain phase-specific output conditions fr< 
Chases I, II, III, and IV, respectivelv. The 



. JLJ.lt: VCU- JLdUJLt: JLU.tJUUJLL.y J-0 

8 2j 
:rom module SQ. 

LUN 19 is used as a "scratch" pad and stores the values g^. of 
lodule SQ for each phase, with each successive phase-specif ic~output 
rector written over the previous vector. 

Information is read from and/or written on the LUNs in the followd 
>rder as a simulation run proceeds (see fig. 3) : 

1) Module Sp: 

Z Functions: ZCOMPiZ/ calls special function TABLES 1. 
Subroutine TABLE31 

This subroutine writes most of table SlL' on LUN 17. The paramete 
hat are local to the lower module Sj_ are written from that module. 

G Functions : LUN 1 has been constructed before the simulation anc 
contains the initial conditions for phase I. Only g^ returns a value * 
:his time, the other g functions produce no (zero) outputs at k 1 =0. 

;2) Module Sp (phase I) : 

Z Functions: ZCOMP, on the initial entry only, makes a call to tt 
ipecial function TABLE32. 

Subroutine TABLE 3 2 

This subroutine writes the portion of tab'ie 3 that contains the 
iestruction/consumption ratios and the larval growth rates on LUN 17. 

F Functions: Function 4 4 makes a call to special function $2- 
Subroutine S2(b) 

This subroutine writes the value of each state variable for each 
ccasion (k 1 ) on a random access file, LUN 25, in binary and a duplicat 
utput on LUN 21 in BCD. The material on LUN 25 is then available for 
.ccess by an external program TABLE4 that will produce detailed infor- 
lation on any state variable at any occasion for each tree class during 
L simulation; this is produced as output table 4 on LUN 20. 



/ Two special-purpose standard FLEX subroutines, ZCOMP and YCOMP, are 
.ssociated with each module in a model. ZCOMP is processed first and YCOMP 
.ast in each time step. ZCOMP handles all input-variable processing, and 
'COMP handles all output-variable processing. 

1 o / 

To distinguish computer output tables from the tables in this paper , 
.he computer output tables will be written in italics, table. Subroutines ai 
;iven in all cap., e.g., TABLE32. 



Y Functions: YCOMP.19/ writes the values of y-^, y 2 , and 73 on LU 
rhese are the year-end values of the variables, that is the phase- 
specific population-branch model outputs. 

Module Si is processed once for each tree class with no calculat: 
}f any g function values in module SQ during this processing. 

(3) Module S n : 

G Functions; Function g 2 , through special function REDIST, writi 
the values g 2 . on LUN 3. 

Function go writes the values g^. on LUN 19. 

Function g^, utilizing the information on LUNs 1, 2, 3, and 19, 
writes new values on LUN 4, where x^ is replaced by the tree class 
specific values of y-j_ from LUN 2, xg is replaced by 3. from LUN 19 ai 
xy is replaced by the values g^. produced by this function. LUN 4 no^ 
contains the initial conditions for phase II. The replacement done b; 
this function constitutues the updating of variables Xc, xg, and xy; 
this is normally an automatic part of the FLEX algorithm, ~but because 
3f the vector manipulations occurring in this module, the updating mu 
De done by a user-defined function (see Vector Processing in FLEX2). 

All other values are simply transferred unchanged from LUN 1 to 
LUN 4. 

(4) Module S 1 (phase II): 

F Functions: Function f^ 4 makes a call to special function S2. 
Y Functions: YCOMP writes the values of y 1$ y 2 , and y3 on LUN 5 

Module S-^ is processed once for each tree class with no calculat 
of g function values in SQ. 

(5) Module S n : 

G Functions : Function g 2 writes the values g 2 . on LUN 6. 
Function g3 writes the values 3. on LUN 19. 

Parallel to what is done in (3) above, function g^ writes the 
/alues g^. on LUN 7, as well as transferring other appropriate inform, 
:ion from LUNs 4, 5, 6, and 19 to LUN 7. LUN 7 contains phase-III 
Initial conditions. 



i q / 

' Refer to footnote 17 on previous page. 



F Functions : Function f^ 4 makes a call to special function S2 
Y Functions : YCOMP writes the values of y-^, 72, and y^ on LUN 

Module S;L is processed once for each tree class with no calcula 
of g function values in SQ. 

(7) Module S Q : 

G Functions ; Function g2 writes the values 2. on LUN 9. 
Function g3 writes the values g~ on LUN 19. 

Function g^ writes the values g^. on LUN 10, as well as trans- 
ferring other appropriate information from LUNs 7, 8, 9, and 19 to L 
LUN 10 contains phase- IV initial conditions. 

(8) Module S-j (phase IV): 

F Functions : Function f^ 4 makes a call to special function S2 
Y Functions: YCOMP writes the values of yj_, y2* and y3 on LUN 

Module Sj_ is processed once for each tree class with no calcula 
of g function values in SQ. 

(9) Module S Q : 

Y Functions; YCOMP calls REDIST for computation of the mean fa 
egg density over the stand after phase IV. REDIST then computes the 
redistributed egg density for each tree class and writes them as Xy 
LUN 12. 

YCOMP then calls the special function subroutines TABLE33 and 
TABLE12 to complete table 3 and produce the output tables 1 and 2. 

Subroutine TABLE33 

This subroutine completes table. 3 on LUN 17 by writing the fall 
mean egg density per model branch over the stand for each phase. 

Subroutine TABLE12 

This subroutine produces the annual input-output summary of the 
model branch for each phase and tree class as table 1 on LUN 15. It 
also determines the postoutbreak probabilities of primary and second 
mortality, probabilities of degrees of top-kill, height- and diamete 
growth reduction factors for each tree class, all according to the 
classification structure provided by Boyd Wickman (1963, 1978a, 1978 



.nd 14. This information is used in producing outbreak impacts as tal 
>n LUN 16. 

Additionally, table 1 is output on LUN 18 in a format suitable fc 
'iewing from an interactive terminal screen. This is done for table 1 
.s well, but with the requirement of outputting it in three parts on 
,UNs 22, 23, and 24. 



STAND-POSTOUTBREAK MODEL 

The stand-postoutbreak model is the filter that takes the defoli- 
ition history from an outbreak simulation and associates with defoli- 
ttion for each tree class the expectations of: direct and secondary 
lortality; degrees of top-kill; and radial- and height-growth impacts. 

The development of the mortality and top-kill expectations are 
>ased on the history of host tree plots associated with the DFTM out- 
reak of 1972-1974 in eastern Oregon (Wickman 1978b). The model branc 
lefoliation is translated into crown defoliation via the branch-tree 
elationship. 

For each tree class, branch defoliation is converted into crown 
iefoliation using the logistic equation 



y = 



p i 



P 2 + P 3 e 



PI X 



rtiere for x,ye[0.1]: 



x: model branch defoliation as a proportion 
of nominal. 

y: proportion of crown totally defoliated. 
P 1 = -22.579678 
P 2 = -22.996253 
P 3 = -598,181,961. 

"hese parameters were obtained by nonlinear least squares curve fit tc 
leans of data from the six crown-defoliation classes of Wickman (19781 
>rown defoliation determines one of seven possible effects vectors. r . 
rector is used to partition the trees of the tree class into direct ai 
secondary mortality or a range of top-kill and growth loss. The set c 
r ectors that determine these losses is derived from the input files 
:ead from LUNs 13 and 14 given in appendix C. The same algorithm is i 



as follows: the first seven contain direct mortality and top-kill 
Information; the remaining two contain secondary mortality information. 
Letting A(i,j) be the jth entry from the ith row, one has the following 
:orrespondences : For i = 

1: all parameters correspond to trees (a tree class) with 
less than 15% of the crown totally defoliated, 

2: all parameters correspond to trees with between 15% and 
of the crown totally defoliated, 

3: parameters for trees with between 35% and 65% of the cro 
totally defoliated, 

4: parameters for trees with between 65% and 85% of the cro 
totally defoliated, 

5: parameters for trees with between 85% and 95% of the cro 
totally defoliated, 

6: parameters for trees with between 95% and 99.5% of the 
crown totally defoliated, 

7: parameters for trees with greater than 99.5% of the 
crown totally defoliated, 

md for each of the above values of i, the following definitions are 
*iven of A(i,j) for j = 

1: expected proportion of trees for which mortality is 
directly attributable to defoliation, 

2: expected proportion of trees with no top-kill damage, 

3: expected proportion of trees with leader-kill, 1 year's 
growth only, 

4: expected proportion of trees with top-kill less than 10% 

5: expected proportion of trees with top-kill less than 25% 

6: expected proportion of trees with top-kill less than 50% 

7: expected proportion of trees with some top-kill, 

8: expected loss in diameter growth from defoliation. 

The last two records in these files contain the information neces- 
ary to remove secondary mortality from the arrays described above, 
'hey are as follows: A(i,j) for i= 

8: expectations of secondary mortality associated with 
bark beetle attack, 

9: expectations of secondary mortality from all other 

causes, e.g., windthrow, winterkill, natural attrition, 



. r or j = 

1: expectation for trees without top-kill, 

2: expectation for trees with leader-kill only, 

3: expectation for trees with more than leader-kill, 
but less than 10% top-kill, 

4: expectation for trees with top-kill between 10% and 1 
5 : expectation for trees with top-kill between 25% and f 
6: expectation for trees with top-kill greater than 50%. 

Because the top-kill information on the external files is cumulative 
: erences are taken between successive classes to get the number for 
L interval. From the resulting classes, secondary mortality is 
ived after combining the two secondary mortality factors that caused 
>ark beetles and that caused by other factors such as windtnrow, 
r , fire, and other "unknown" causes. 

The two forms of secondary mortality are assumed to act indepen- 
:ly, hence the two vectors are combined by the rule: 

(m- . + m . - m 1 .m ,) =v.. 
1,J 2,j l,j 2,j' j 

scalar- or dot-product of this vector with the input vector of 

.duals of direct mortality is then used to compute secondary mortality 

each defoliation class: 



i r^ * is the vector of residuals (ith record of LUN 13 or LUN 14) , 
si is the derived secondary mortality that will be used in table 2 
simulation output (appendix C) . 

After secondary mortality has been subtracted from the residuals of 
ict mortality, the remaining numbers are the expected proportions 
: are to receive various levels of top-kill, and one class that 
.cates the proportion receiving only growth reductions. 

Growth reductions attributable to defoliation are changes in 
leter growth at breast height and height growth. Growth reduction is 
inction of maximum model branch defoliation. These assumptions 
tit in the following table of possible effect vectors. One vector is 
len by species and defoliation. 

An example of defoliation impact summaries is given in table 2 of 
output section of appendix C. 



Ou tput 



of 



** 



* 



uglas-fir host: 



Expected top-kill 



Percent 
of crown Direct 

totally mortality 
iefoliated 



0-15 
16-35 
36-65 
66-85 
86-95 
96-99.5 
99.5-100 






.0090 
.0280 
.1730 
.4770 
.9230 



Secondary 
mortality 



.1127 
.1229 
.0927 
.1000 
.1950 
.1070 
.0158 



Top-kill 
(total) 



.1129 
.2081 
.4032 
.4820 
.2929 
.1864 
.0274 



Growth 

impacts 

only 



.7744 
.6690 
.4951 
.3900 
.3391 
.2296 
.0338 



Leader 
kill 


1-10% 
of 
crown 


10-25% 
of 
crown 


25-50% 
of 
crown 


50-90% 
of 
crown 



.0675 
.1360 
.1744 
.2376 
.1400 
.1422 



.0454 
.0696 
.2058 
.1606 
.0607 
.0359 



.0209 .0053 
















.0230 

.0838 

.0782 









.0026 



.014C 
.008: 
.0012 



Grand fir host: 



Percent 

of crown Direct 
totally mortality 
defoliated 



Expected top-kill 



0-15 
16-35 
36-65 
66-85 
86-95 
96-99.5 
99.5-100 






.0090 
.0280 
.1730 
.4770 
.9230 



Secondary 
mortality 



.0451 
.0422 
.0455 
.0592 
.0957 
.0790 
.0116 



Growth , A I-* 10 - 25% 25 - 50% 5 ~ 9C 
Top-kill ? s Leader f of O f of 

(total) f* kill crown crown crown 



.0544 
.1261 
.2674 
.3817 
.3387 
.2137 
.0315 



.9005 
.8316 
.6781 
.5311 
.3926 
.2302 
.0339 



of 
crowr 



.. 

.0331 


.0123 


.0039 


.0029 


.002; 


.0877 


.0294 


.0068 


.0019 


.000- 


.1322 


.0882 


.0317 


.0113 


.003' 


.1723 


.1022 


.0424 


.0480 


.016 


.1893 


.0580 


.0217 


.0157 


.054 


.0963 


.0401 


.0274 





.049 


.0142 


.0059 


.0040 





.007 
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APPENDIX A VERSION 3.1 SOURCE CODE 

Page 
Command file 57 

DFTM stand-outbreak model 

Module S Q 58 

Module S- 61 

Stand-postoutbreak model 

Module S 65 

Module S 73 

Subroutine TABLE12 flowchart 74 



00 ON INPUT UNIT 
3CCG1 -MQQUL-GETO=1 
00002 -NUM3ER=7 Q = U 

OCG03 -TMAX=I FUNLOA o=*\/3iS3 

OOG1H -YMAX-3 NODUM^ NOLOG TTYQFF NOSUMMARY 

OC005 -TITLE= 5(1) R\/ISFO STANT OUTBREAK MODULE 

00006 -3 O=0. , 3* , &., 0., w 0. , _C, ^, f C . ^ C., G., 

o oc 0-9- - o^ f 5^4 $ 1 , -b* 1 3 * o b i l^iSo ilololi^aloozlcloni^oloio* 

JOC1G -0.016,0. 3*2,0 . uC5,G. CO 6,0.007, Q. 021, 0. 033, 0* 056,G .5, 0.62, 
00011 -C --------- 



G0013 -RO= l.J, 0.^5, 0.1, C., C., G.t 1 . L. , 0.6, 0.5, G,l< 
OCr-1'4 - C., 3., D.G2, 3.02, O.G2, G.C2, 0.02, 0.02 



OC015 - 1 

00016 -NUMPER=* 0=10 TAP.GF T 

OCC17 -TMAX-1 FUMLOAO='*V31S1 

00^1* - U M D H J N = ^ 

0001Q -0 = YU) C=Y(2) G = Y(^5 

00028 YMAX=3 

00021 -LP=Y(1) LP=Y(2) L D =Y(3) 

OCG2? -0=Y(1) C = y (2) D=Y(3* 

C r! 23 ""BCI *32, . J ? , e^c, . ^ 2 , OP, .02, 5.^-, 6.25, 6.2 F, 2 

0002^ - 2.27, 2.5, 5.^4, ^>?^5, 6.25, 3.69, 3.29, 3. , "l. 19 , \ . C 8 

OOC2F -TlfLE = siMULAT ICN ON*>HE RAISED MOOEL^VFRS . ^3 / " 
OCG2 7 -T T YQFP NOSUM^APY 



OOC29 -FIRSTSU5 O = l,l 

000 30 -MO n UL EGET( ) =1 9 1 
OUt 31 -GHOST C) =1 
TOT32 -INPUT=6C 
EOF CM INPUT UMT 



[JB3 -C 



005 - 


{ 


006 - 


V 


00? - 2QFQ ? 


008 - 


t 


009 - 


\ 


01 -C 




Oil -C 


* * \ 


012 -C 




013 - 


\ 


014 - 


1 


015 - 


] 


016 -C 




017 -C 




018 -C 




019 -C 




02C -C 




1021 - 




022 - 




1023 - 




024 - 




1025 - 




1026 -C 




1027 -C 




1025 - 




1029 - 




1030 - 




3031 - 




3032 - 




3033 - 




103^ - 




3035 - 




)036 - 




3037 - 




3038 -C 




3039 -C 




3040 - 




0041 - 




0042 - 




0043 - 




0045 - 




0046 - 




OH47 - 




0048 - 




0049 - 




0050 - 




0051 * 




0052 - 




0053 - 




005^ - 




0055 - 




0056 - 


1000 


0057 - 




005* - 


1001 


0059 - 




0060 - 




0061 - 


1002 


{Jfi6 2 


10 


0063 - 




0064 - 




0065 - 




0066 -C 




0066 " 




IQT169 - 




I007C - 




10071 - 




10072 - 




10073 - 




0074 - 




10075 - 




10076 - 




10077 - 




10079 - 




10079 - 




10080 - 





SUBROUTINE ZCOHP fK ,X jB.R.Zj TO , Tj G 

DIMENSION XCllfBCllfRCl),ZCl)iGlll 



WRITEC61,2000) 
FORMAT(# ENTFR 



NUMBER OF TREE CLASSES IN STND#) 



E THE BEGINNING OF TABLE 3 *** 



CALL TABLF31(e,R) 

RFTURN 

END 



FUNCTION GOKKX,B,R* ZfKP) 
niHFNSTON XU* ""* " 



RETURN 
END 



FUNCTION G2 CK, 
DIMENSION XU) , 

IF 2 7i<P.ECUlPHASE) RETURN 
IF CKP.EO.Oy RETURN 



,R 



CALL REOISTtKPtINUM^X) 

IPHASF=KP 

RETURN 

FND 



FUNCTION r,03(K,X ? B ? R,Z 
DIMENSION X(l)tB<l)tR< 



(l) 



tZ(l) 



TO - 

ir"<TP.EQ.IPHASE) RETURN 
IFCKP*EQ* 0) RFTURN 
NN=IFIX(7(D) 



LUN1=LUN*1 



5) 
X6 



REWIND LUN1 

REWIND LUN2 

REWIND LUN19 

DO ID J=1,NN 

READ(LUNlflOOO) 

FORMAT<24X* 2F11 

RFAD(LUN?tl001) 

FORMAT(1?X,1F12.6) 

Gi=FMINC<l-Z?/100)*Z3,X6> 

WRITECLUN19,1002) G3 

FORMATC1F12.6) 

CONTINUE 

IPHA.SE=IP 

RETURN 

END 



FUNCTION G04(K,X,R, : R,7,KP) 
DIMENSION X<1),3(1) ,*M1) ,Z(1) 



IP=KP+1 

IF (IPECUTPHASE) RETURN 
IF (KP^EQ.O) RETURN 
NN=IFTX(Z(1) ) 



LUN1-LUN*! 
LUN2=LUN*2 
LUN3=LUN-f3 



LUN19-19 



00062 
00083 

OQ08 

00086 

00087 

ooose 

QQ089 
00090 
00091 
00092 
00093 

00095 

00096 
00097 
00096 
00099 
0010C 
00101 
00102 
00103 
0010** 
00105 

00107 

ooioe 

0010? 
00110 
00111 
00112 
00113 

00115 
00116 
00117 

oone 

00120 
00121 
00122 
00123 

H0125 
00126 
00127 
00128 
00129 
00130 
00131 
00132 
00133 
0013<* 
00135 
00136 
00137 
0013P 
0013C 



RFWINO LUN? 
REWIND LUN3 
REWIND LUN** 
REWIND LUN19 
00 120 I=1,NN 
*FAO (LUN1 2 2000) LL 

- 2000 FORMAT <U2 f UF11.5t 

READ (LUN2,22QO) X5 
READ CLUN3 f 2200) G2 
READ (LUN19,2?OOJ G3 

- 2200 FORMAT (1F12*) 



,113) 



X7=(1.-B<IPP) )*G2 

- 120 WRITE(LUN**,2SQP) IZ6, Z*, Z5, Z 2, Z3,X5 ,X6 t X7 W 

- ?300 FORMAT 1 1 12, 7F11. 5, 113) 

RFWINO LUN4 
IPHASE=IP 
RETURN 
END 



-C 
-C 



-C 
-C 

-C 

-C 
-C 

-C 



-C 



SUBROUTINE YCOMMK* X B ,R Y f XS,T) 
DIMENSION X(1J fBCll ,R(l) f Y(l) f XS(1) 

IR57=IFIX (8(57)) 

CALL PFOISTCKP,IB57,B,X ) 

*** FINISH TABLE 3, ON LUN 17 

CALL TABLE33CX,B) 
*** WRITE TARLFS 1 ANO 2 *** 

CALL TA8l12(P) 

PFTURN 

FNH 



-C 
-C 

-C 



-C 

-C 
-C 

-C 
-C 
-C 



001<*1 - 



001^5 - 



SUBROUTINE REOIST(KP,INU^,eX) 
DIMENSION 8(11 ,Xf 1),XYCOMPf l 

INITIALIZE LUN VALUES ANO REWIND 

LUN1=(KP-1)*3-1 
LUN2=(K D -1) 
LUN3=CKP-j) 
RFWIND LUN1 
RFWIND LUN? 

^ REDISTRIBUTION OF FGGS ** 

INITIALIZE LOCAL VARIABLES" 

SUM=TQTSTFH=00 

CALCULATE TOTAL EGGS OVFR STAND (SUM) 

AND TOTAL WEIGHTfO SUM OF TREES I TCTSTEMI 

00 111 J = 1,INUM 
RFADCLUN1 ,1001) Z^,7^ 

RFAD CLUN^lOoS)' X7 
1000 FORMAT (2<*X,F12.6) 



- 1<] 

-C 
-C 
-C 
-C 



-C 
-C 



CONTINUE 



CALCULATE X(KP) WHERE KP IS PREVIOUS PHASF 
STORE LOCALLY IN XYCONP(KP) UNTIL PHASE IV 
TO PPEVFNT XS=n FROM ERASING STORED VALUES 



XYCOMp(KP) =S 
REWIND LUN2 



CALCULATE VALUES OF G<2,J) FOP SO HERE, 

INSTEAD OF IN FUNG 1ION GO? CODE AND WRITE ON LUN 



DO 20 J=1,INUM 



) - PEAO CLUN?, ioooj XT 

. - G2 = X7+B(62) MXYCOMP(KP) -X7) 

! - WPITE <LUN3 f lGll3J G? 

! - 1010 FORMATCF1Z.6) 

- 20 CONTINUE 

J - IF(KPEQff) GO TO 26 

: - RETURN 

' ~C 

! -G TRANSFER VALUES FROf LOCAL STORAGE TO X ARRAY 

I -C 

\ - ^0 CONTINUE 

! - RETURN 

:N INPUT UNIT 



uouu* -u 


l-LtA i-UNVlIUNb t-UK V5l Ur \ W NULJULt bllJ 


00003 -C 




0000* -C 




00005 - 


SUBROUTINE ZCOMP CK* X,3* R,Z* TO tT t G) 


00006 - 


DIMENSION Xm iB(lltR(l),Zl),G(l) 


00007 -C 




00006 -C* 


** READ IN THE Z VECTOR, INITIAL X VECTOR *** 


00009 -C 




00010 - 


CALL VARIN f?(3),G(D) 


00011 - 


IP=IFIX <Z(3)) 


00012 - 


LUNl = fTP 1) * '3 + 1 


00013 * 


READ CLUN1, 1000) ISPECTESt Z f 1) ,Z(?)Xfl),X<2) f X(3) 


0001*4 - 1 


000 FORMAT ( 112, 22X , 5F11 . 5) 


OOQ1E - 


Z C t+ ) =Fl. OAT ( IS PEC IES) 


00016 - 


X f f+ ) =0* 


0017 -C 




oooie -c 


*** SET PHASE-SPECIFIC PARAMETERS B(2o> AND 3(3436) 


0002C - 


B (29) -7(^) 


00021 - 


CALL VARIN(B(3lf8tIP*i8)) 


00022 - 


CALL VARIN (B ( 35) , B I IP+521 > 


00023 - 


CALL VARIN IB (36 1 ^ B (IP+62) ) 


02^ *C 




00025 -C* 


** REAH IN HOST AND PHASE SPECIFIC MORTALITY VECTORS 


00026 -C 




00027 - 


DO 6 J=l,6 


00028 - 


I1=J+(IP-1) *6 


00029 - 


T2 = J + (IP *3) *6 


00030 - 


j^ = j4-(jSP EC IS 1 ) *6 


00031 - 


CALL VARIN < P ( J) ,R( J4-12) ) 


00032 - 


CALL VARTN <R(J)*BCI1)) 


00033 - 


CALL VARIN (R ( J+6) , R ( 12) ) 


QOQ34 - 


6 CALL VARIN ( R ( J+ 1?) , R (I 3) ) 


00t335 - 


IF (K.NE.O) RFTURN 


00036 -C 




130037 -C 


*** FINISH TABLE 3 ON LUN 17 *** 






0003? - 


CALL TABLE32C8) 


0*4-0 


RFTURN 


*f 1 ~ 


END 


000^2 -C 




Q00*f3 -C 




Q0*5 -C 




o .0 n ** - c 




000^7 - 


FUNCTION G01 K^ X, B 9 R t Z KP) 


000*4? - 


DIMENSION XC1) ,Rfl) tRd)fZ(l) 


000^9 - 


G01-0. 


00050 - 


IFCKP.Ed* B) G01-(Z?)*Z(l)/100.)*{CXfll-t-X<2)*7C2))/ 


00051 - 


A (? : 0*Z(2))~G'(19)*^(27M 


00052 - 


RFTURN 


00052 - 


FND 


00054 -C 




0005* -C 




OOH56 - 


FUNCTION G03 (K, X, B 9 R * Z* KP) 


00057 - 


DIMENSION Xll1B(l),R(l)tZCl) 


00058 - 


G 3 - 


00059 - 


IF' (KPEQ.8) G03=X(1) 


00060 - 


PFTURN 


00061 - 


FND 


00062 -C 




00062 -C 




0006* - 


FUNCTION G0 ( K, X B fR,Z,KP) 


QOQ6* - 


DIMENSION X (11 , 8(1) ,R(1) Z(1) 


00066 - 


G 1+ ~~ 


00067 - 


IF (KP-GT.O.ANP.KP.LE.6)(;OJti-FMINCGfl6) ,XCD) 


00068 - 


RFTUPN 


00069 - 


END 


0007C -C 




00071 -C 




00072 - 


FUNCTION G05 ( K , X 0^ R iZ KF) 


00073 - 


DIMENSION X(l) ,B(l)f R(l) f ZC1) 


00075 - 


TF 5 TKP?GE.1AND.KP.LE.61 G 05 = FMIMG (16) tX (2) ) 


00076 - 


RETURN 


00077 - 


END 


00078 -C 




OOH79 -C 





-c 
-c 



-c 

-c 



-c 
-c 



-c 



-c 
-c 



-c 
-c 



=. 
IF<KP.GT.(UANO.KP-LE6) 606= (l.-R (KP) )*< 1 ,-R <KP*6) )*<!.-* 

RETURN 
END 

FUNCTION Gfl7 ( K f X, 3,R , Z? KP ) 
DIMENSION Xfl)fB(l)R<l)f7Cl) 

IFfKPlGT.QANOKP.LE,61 G 07 = 1. -RCKP + 12 ) 

RETURN 

END 

FUNCTION G08 ( K, X, Bt *i Z KP) 
DIMENSION X(l) 8tl) t*<l)*7<l) 

IF 6 7KP*EQ.O.OP.KP.GT.6) RETURN 



RETURN 
END 



FUNCTION G09 (K,X, 8 ,?f Zt KP) 
DIMENSION X(l) ,B<11 tR(l) tZCl) 

P.EQ,n.OR.KP,EQ^) GTOT-=0 



G09=GTOT 

RFTURN 

END 



FUNCTION GtO ( K, X* B,R f Zt KP) 
DIMENSION X<1) fB(l) tR(l) tZCl) 



- 10 RETURN 



IFCKP.EQ.' 
IFtKP.EQ. ' 
IFCKP.EQ*' 
RETURN 
END 



FUNCTION Gil CK,X,8fR* ZtKFI 
HTMENSION X ( 1) B (1) ,* (1 ) Z ( 1) 
Gll=(Gf 8)*G(10))*X(^) 
PFTURN 

END 



FUNCTION G12 C K, X f 8 , R f Z f KP) 

DIMENSION X(l),B<l),Rtl),7(l) 

G12=0. 

IF(X(3) .1E.Q..ORX(1) +X<2) .LE.O.) GO TC in 

IFCKP.GE.l.AND.KP.LT.U) ~ "* ^""" % 

IF (KP.EQ.^) G12=P(23) 

IFCKP.GT.if.AND.KP.LE.B) 



-c 
-c 



FUNCTION G13 C K f X B,R f Z f KP] 
DIMENSION XC1) fB(l) tRCDvZf 1) 



IFCKP.EQ. 0) 



IF(KP*GT.O.AND.KP.LE.6) G13-X ( f I *(EXP( 1 0. *G t 12) *PHI) 

RFTURN 

END 



-!) 



FUNCTION Gl* (K t X^8tRZt 
DIMENSION X(l) t BCl) ,R(1 )tZ(l) 



RETURN 

FNO 



-c 
-c 



00162 
00163 
0016* 
00165 
00166 
00167 
0016? 
00169 
00170 
00171 
00172 
00173 
0017** 
00175 
00176 
00177 
00178 
00179 
00180 
00181 
001*2 
00182 
0018* 
00185 
00186 
00187 
00186 
00189 
00190 
00191 
00192 
00193 
0-019* 
00195 
00196 
00197 
00198 
00199 
0020C 
00201 
00202 
00203 
0020* 
00205 
00206 
00207 
00208 
00209 
00210 
00211 
00212 
00213 
0021* 
00215 
00216 
OQ?17 
00219 
00219 
00220 
00221 
00222 
00223 
0022* 
00225 
00226 
00227 
3022f 
00229 
QO?3H 
00231 
0023? 
00233 
0023*4 
00?35 
00236 
00?37 - 

00238 - 

00239 - 
002*0 - 
002*1 -C 
002*2 -C 



-C 
-C 



FUNCTION G15 (K, X, 8 ,R , Z ,KP) 
DIMENSION X(11 ,8(1) ,R(1),Z(1) 

IF(KP.EGU O.OR.KP.GT.6) GO TO 15 

DO 10 1=0,10 

ALPHA=FLOATI1Q-I) 

G15=ALPHA 

TF(UVl(KP,X f 8,ALPHA) .LT .X(D) RETURN 

CONTINUE 

RFTUPN 

G15-0. 

RFTURN 

END 



FUNCTION G16 <K ,X t B ,R , Z,KP) 
DIMENSION X <11,B(1) fR(l)f?(l) 



, 

IFCKP.EQ* 0OR.KP*GT.6)RETURN 
/\LPHA=PMIN(G( 15) +1* , 10.) 
IFCX fl) *GT. 0. )G16=U VI CKP,XtBt ALPHA) 
RETURN 
END 



FUNCTION G17 ( K, X,B,R , Z, KP) 
DIMENSION Xm,B<l)*im),Zfl) 
IF(KP.ECUO.OR.KP.GT.6) GO TO 20 
G17=10.-G(15I 
ALPHA=G(15) 



DO 10 T=0 ,N 
ETA=FLOAT (II 
ETAPLUS=ETA-f-l. 

IF( < UV2 ( KP^X * B, ALPHA, FT A PLUS) -fX11-UVl(KPfX,B, ALPHA) ) ) 
r.GT*X(2)) GO TO 15 
CONTINUE 
RETURN 
G17=ETA 
RFTURN 
G17-0* 
RETURN 
END 

FUNCTION G18 I Kt X, B ,P , Z ,KP) 
DIMENSION X(l) ,8(1) ,R<1) ,Z(1) 

G18=0. 

IF (KP. EQ.Q*OR.KP.GT.f ) RETURN 

ALPHA=G (15) 

ETA=10-ALPHA 

G18=UV? (KP,X^e, ALPHA f FTA) - tXtll-UVK KP,X,B*ALPHAM 

IF (Glfl.LT. 0. ) Gia=fU 

PETURN 

FMO 

FUNCTION G19 ( K, X, B,R , Z, KP) 
DIMENSION X(l) ,BC1) ,RC11,7(1) 

IF(KP.FQ.O*OR,KP.FQ .9) CTOT=0. 

IF(KP.GI. 0. ANO.KP.LE.6)GTOT=GTOT*(10.- < G ( 15 ) +G ( 17) I) 

G19=GTOT 

RETURN 

END 



FUNCTION F01Q1<K,X,R,R,Z,KP) 
DIMENSION X(l) ,8(1) ,Q(1) ,Z(1) 



RETURN 



FUNCTION FQ10?(K,X,B,R,Z,KP) 
DIMENSION Xtl),3(l),3(l) ,ZC1) 



RETURN 

END 



FUNCTION Ft!2Q2CK,Xf P,R,Z,KP) 
DIMENSION XC1) *Bdl **d),Z(i 



RFTURN 
END 



FUNCTION FOSOSiK^XiBiRiZiK*) 

DIMENSION X<1* ,8(11 tRUltZCl) 

F0303=in. O'Gdll 

RETURN 

END 



FUNCTION FIHOt* ( K f X f 8 f R,Z KP) 

DIMENSION XC1) f9fi) v 

FO<0<*=G(13) 

CALL S?(X t 

RETURN 

END 



SUBROUTINE YCOVP ( K, X , B, R 5 Y XS f T) 
DIMENSION Xfl) ,^(11 9 P(l)fYCl) 
D 1 ^ 10 T = l,3 
10 Yd) =X(I) 

IF (K.EQ.O) RETURN 
IP=IFIX(B(29)) 



WRITP (LUN2tl0001 Y ( 1 ) * Y f 2) , Y ( 3) 
- 1000 FORMAT (3F12.6) 
RFTURN 

END 
-C 
-C 

FUNCTION UV1 (KPtX,8, ALPHA) 

DIMENSION XC1) ,Bd) 



IFCKP.EQ. O.OR.KP.GT.6) RETURN 



= U*(ad9l+E(20)/Gd2))MEXP<Gd?>)-l.)*Q.001/D(21>*B<!>* 
C<1*- (G(6>*XP(G(12) ) ) **ALPHA)/(1.-G (6) *EX^ CG d 2) ) ) 
RETURN 
FNO 
-C 
-C 

FUNCTION UV? <KP,X,8,ALPHA,ETA) 
DIMENSION Xd) ,Rd) 



IF(KP.EQ.O.OR.KP-GT.ft) RETURN 



IFCG(T) *EQ.O. .flND.ET/UEQ. C. ) RETURN 

UV2 = U*(B<191 +fl(20)/G(i2M(FXPIGd?M-l.)*0001/8(211*Bm 
CG (6) **ALPHA*EXP(ALPHA*Gd21 ) *< 1 .- ( G ( 6) *G C 7 ) *FXP ( G (1 2) ) ) **ET A) 
Cd.-G (6)*Gf 71 

PETUPN 

END 
I INPUT UNIT 



00002 - DIMENSION 8(11 ,R(1) 

00003 -C 

0000** -C BFGTN TABLE 3 ON LUN 17 

00005 -C OBTAIN T I ^F AND DATE 

00006 -C 

00007 - TIME^G* 

00008 - DATF-=0. 

00009 - EM = Q 

flOQlO - CALL XREQ ( 339 , , 0, EM , T IM ) 

00011 - CALL XREQ ( 32 , , 0, EM, DATE) 

00012 -C 

00013 -C WRITF HEADING ON LUN 17 
0001^ -C 

00015 - WRITE (17,1000) TIMF,DATE 

00016 -C 

00017 - 1000 FORMAT(*l*,/*-*,/*~#,9X ,*TABLE 3? PARAMETERS FOR 
OOOie - C,* OUTBREAK-'*, 

00019 - C* MORTALITY, FOLIAGF DESTRUCTION, AND OTHER.* 

00021 - C*HOST TREE, STAND* PHASE, OR INSTAR SPECIFIC*; 

00023 - C*PUPAE*,5X,*OVERWINTFR*,6X,iEGG^MASS*/* * ,,-.,. 

0002^ - C6X,*3*, 6X,****, 6X, *5*, 6X,*6*, 9X, MORTALITY MORTAL] 

0002 C - H 9X,*SIZE*/* *, 

00026 - C135(*-*)/*0 GAILY MORTALITY RATE*/* *,2M*-*)) 

00027 -C 

00028 -C WRITE BACKGROUND DAILY MORTALITY RAT r S 

00029 -C 

0003C - WRITF (17,1010) (R( J) , J=13 f Ifi) 

30031 - lOl'O FORMAT(*0*,9X, *BACKGR OUND*, 8X, *I, I T , II I, I V OF,GF*, 

00032 - C6(1X,F5.3 ,1X1) 

00033 -C 

00t?3^ -C WRITE PATLY DISEASE MORTALITY RATES 

0003*5 -C PHASFS I, IT, III, IV 

00036 -C 

00037 - WRITE(17,1020) (R(J),J=l,e) 

00039 - " WRITE ( 17 , 1 30 ) (8 ( J ) , J^7 , 1 2 ) ' >,-,, 1 
---- - - 1Q30 FORMAT(* * , ^OX , *II *, ?X , *DF ,GF* , 5X, 6 ( IX , F5 . 3, IX ) ) 



FQRMAT(* *,30X,*ITI*,6X,*OF,GF* t 5X,6(lX,F5. 3, IX) ) 

000 1* 3 - WRITE (17,10^5) ( P ( J) , J =19, 2<t) 

000^ -10^5 FORMAT (* * , 3 X, *I V* , 7X , *DF ,GF* ,5X , 6 (1 X, F5 . 3, IX) ) 

OtSO^ii -C WRITC" DAILY PAR ASTTE/PREDATOR MORTALITY RATES 

000^7 -C CHASES I, II, III, IV 

ooo^e -c 

OOPU9 - WRITE (17,1050) ( B( J) , J-25 , 3 0) 

00050 - 1050 FORMAT(*0*,9X ,*PAPA SITE/PREDATOR* ,5X, *I#, 7X , *OF, < 

00051 1 6 (IX, 1F5.3,1X) ) 

0005? - WRITFd 7,1060) ( B( J) , J=31, 36 ) 

00053 - 1060 FORMAT(* *, 30 X ,*II* , 7X , *DF , GF* , 5X , 6 ( IX ,F5 . 3 

0005^ - WRITE (17,1070) ( 9 ( J ) , J = 37 , 42) 

00055- 1070 FORMAT(* * , 3fl X , * II I * , 6X , * OF , GF * ,5X, 6 (1 X , F5 * 3, 1 X) ) 

00056 - WRITE (17,10751 ( 8 ( J) , J = 43 , k 6) 

00057 -1075 FORMAT <* *,3CX, *TV*,7X 1 *DF t GF*,5X,6(lX, 1F5.3, IX)) 

00058 -C 

00059 -C W^IT r HOST-SPECIFIC DAILY STRESS MORTALITY Rfll 
OOH60 -C PHASFS I, II, III, TV 

00061 -C 

00063 - WRITE(17,1Q80) (R(J),J=1,6) 

00063 - 1080 FORMAT(*0*, <5X,*STRESS*, 12X,*I, II, I II ,1 V*, 2X , *DF* ,7) 

0006^ - WRIT^ (17,1090) (R ( J ) , J =7, 12) 

0006^ - 1090 FORMAT(* *, 27X , *I, 1 1 , III , I V t , ?X , *GF * t 7 X , 6 (X , F5 . 1 , X 1 

00066 -C 

00067 -C WRITE SUP-HFOING 
0006^ -C 

0069 - W^ITE ( 17, 1095) 

0007C - 1095 FORMAT(t , 2^ t *- *) / * OCCASION MORTALITY RATE*/* 

^10071 -C 

00072 -C WRTTF PUPATION MORTALITY ^ATES 

00073 -C PHASES I, II, III, IV 
0007^ -C 

OOOT5 - WRITEtt 7, 1100) 8(^9) 

00076 - 1100 FORMAT(*0*,1 OX,*PUPAE*,16X,*I*,7X,*CF,GF*,5 C 7X,F5.3] 

00077 - WRITF (17,1110) R<5(1) 

OOD7P - 1110 FORMAT(* *, 30 X ,* II * ,7X , *DF ,GF* , 55X^F5* ?) 

00079 - WRITEU7, 1120) R(51) 



00080 


- 11?0 


FORMAT** *, 30 X ,#TII*,6X ,*DF,GF * ,55X , F5*35 






00081 





WRITE (17,1125) P(5?) 






00082 


-1125 


FORMAT it #,3 OX, *IV*, 7X ,*OF, GF* , 55X , 1F5* 3) 






00083 


-c 








0008*1 


-c 


W?IT r OVFSWTNTER EGG MORTALITY RATES 






00085 


-c 


AFTER PHASES I, IT, III, IV 






00086 


-c 








00087 


. 


WPITE (17*1130) Rf58) 






QQ08S 


- 1130 


FQRMAT(#Q*,9X, ^OVERWINTERS, 9X,*I - II^,5X,#D 


F,GF/,67X 




00089 





WRITE(17, 11^0) 3(59) 






00090 


~ 11 ** 


FOR^AT( #,?7Xf#II - III#t^X,^DF,GF< f 67X,F5 


3) 




00091 


. 


WRITE (17,1150) R(60) 






00092 


- 1150 


FORMAT(# #,27X,^III - IV OF , GF^ ,F>7X ,F5, 3 ) 






00093 


- 


WRITE (17,1155) 6(61) 






0009*4 


11 c >5 


FORMAT (# ^,27X,#IV- *, 5X , #OF,GF# ,67X , IF 5 


3 ) 




OOH95 


-C 








00096 


-C 


WRITF NOMINAL EGG MASS SIZE 






00097 


-C 


PHASES I, II, III, IV 






00098 


-C 








00099 




WRITEC17, 1160) B(53) 






00100 


- 1160 


FORMAT(# gt,2^C<-#l/#0 NOMINAL FGG MASS SIZ 


F*, ex, 




00101 




")3CI#^7X, #0^, GF #,82X, lrc >l) 






00102 


~ 


WRITE (17,1170) B(54) 






0010 2 


- 1170 








(3010** 




WPITE(17,1180) 3 (55) 






00105 


- 1180 


FORMATC* it*30X ,^ 111^,6 X ,#DFtGF* f 62X, F5 . 11 






00106 




WRITE (17,1190) B(56) 






00107 


- 1 1 9 


FORMAT (^ix tn v t 7X ^DF GF5t 8?x IFS i) 






noios 




RETURN 






00109 


- 


EN9 






FOO CN 


INPUT 


UNIT 






00001 





SUBROUTINE TA8LE33(X,B) 






00002 





DIMENSION X(l) ,B(1) 






00003 


-C 








00004 


-C 


COMPLETE TABLE 1 BY WRITING THE PHAS^-SPEC IFIC ME/ 


00005 


-c 


FGG DENSITIES FOR PHASES I, II, III, IV 


, AND THE 




oouoe 


-c 


ANH THE OUTBREAK RECISTRI BUNION COEFFIC 


IENT 




00007 


-c 








00008 


- 


WPITFU7,3QOQ) (X<I), !=!,**) 






00009 


- 3000 


FORMAT(^0^,135 (*-*}/*- PHASE SPECIFICS, 






TI oo in 


f 


;# MFAN EGG DENSITIES: I=^,F7.3,* TI=*,F7 


* 3 , 




30011 




/^ I II -I , F 7 3 , t IV = ^,1F73) 






00012 





WRITF "(17,29991 B(62) 






00013 


- 2999 


FORMAT (^REDISTRIBUTION COEFFICIENT FOR THF#, 


00014 




:t OUTBREAK? *,F7.3) 






00015 


. 


RFTUPN 






000.16 


- 


FND 






ECO CN 


INPUT 


UNIT 






00001 





SUBROUTINE TA8LE12(B) 






00002 


-r, 








OC003 


-C 


THIS ROUTINE WRITES THE TWO SPECIALLY FOPMATTED 01 


00004 


-n 


SUMMARIES REFERRED TO AS TABLES 1 ANO 2 







00005 


-r: 


TABLE 1 GIVFS THF INPUT/OUTPUT VARIABLE 


VALUES, 


B^ 


noooe 


-C 


OVER THF COURSF OF THE OUTBREAK. TARtF 


? CONDENSE 


00007 


-C 


INFORMATION FU^THEf? AND GIVES THF MORTALITY ANH 


TC 


00008 


-c 


EXPLCTIONS WHICH COMPRISE THE STAND OS 


T-OUTBREAK 


00009 


-c 








0001C 




DIMENSION 8(1) 






00011 


- 


DIMENSION I SPECIES (6) ,IC(15, 3) , A (15, 26) ,C (IE 


,5) , AMORT 


( 1 


00012 





DIMENSION WT(15) 






00013 





DIMENSION IPHASE (4) ,np( 2, 12) 






00015 


. 


DATA (( ISPECIES( T),I=1,6)=^OOUG^,#LAS-#,*FIP 


*,* GRA# 


, 


00316 


r 


* 5fND F/,/IR t) 






00017 


-C 








oooia 


-0 


OBTAIN TI^F AND DATE 






ooni9 


-c 








OOQ2C 




T IMF = o 






00021 


- 


CAT r - n 






00022 


- 


F Mr: 






00023 


- 


CALL XPFO(33B,0,QE^,TIMF.) 







00025 -C 



CALL XREQ(^?PtO 0,EM,OATE) 



88(1, J) 



xjj ,j~i., 
J=7,,.. 



.9r, i-o, .. 

t 12f 1=1,2 CONTAIN 



L NI M i r 
LAST 2 



RECORDS 



i VARIES WITH PERCENT OF CRGWN DEFOLIATED? 

1 OP 8HQ-15Z? 2 OR 91I16-3T/:: 3 OR ID 1 136-65'/.; 

^ OR 11*266-85%; 5 OP 12 * 86-95X; 6 OR 13* t 96-99.5X5 

7 OR 1** t99.5-100X 
VALUES READ IN FOR AMCRT* 

AMORT <Tt n-EXPECTATION OF DIRECT MORTALITY 
FOLLOWING 0R CONDITIONAL EXPECTATIONS, GIVEN 
SURVIVAL OF DIRECT MORTALITY 
AMOPT<I,2)=NO TOP KILL DAMAGE 
AMORT<I,3)=LEADER KILL 
AMORT <I,i*) =1-10X TOP KILL 
AMORT(I*5)=l-?57. TOP KILL 
AMOPT(I t 6)=l-50X TOP KILL 
AMORT(It7)=l-99X TOP KILL 

VALUFS OF AMORT AS RECALCULATED FOR OUTPUT 



2): 
3)' 



AMORT 
AMORT(I, 
AMORT (1, 
AMORT <!,*); 
AMORT (1,5) 
AMORT (1,6 ) 
AMOPT (1,7) 
A TORT (I, 8)- 
AMORT (I t 9): 



:EXPECTATION 
^EXPECTATION 
EXPECTATION 
EXPECTATION 
EXPECTATION 
EXPECTATION 
EXPECTATION 
EXPECTATION 
EXPECTATION 



OF DIRECT MORTALITY FROM DFTM 

OF SECONDARY *ORT., ALL CAUSES 

OF TOP KILL, ALL DEGREES 

OF GROWTH REDUCTION ONLY 

OF LEADER KILL CNLY 

OF 1-10X TOP KILL 

OF 11-252 TOP KILL 

OF 26-507. TOP KILL 

OF 51-99X TOP KILL 



FOR PR ARRAY I=t SIGNIFIES SECONDARY MORTALITY 

DUF TO BARK BFETLE 

1-2 SIGNIFIES S^CONCARY MORTALITY CUE TO OTHER CAUSES 

B8(T,J=1 OR 7)=EXPECT. OF SECOND. MOPT., GIVEN 

NO TOP KILL DAMAGE 

PR(T,J=2 OR 6)=EXPECT. CF SECONC. MOPT., GIVEN 

LFADFR KILL 

eRCT,J=3 OR 9)-EXPFCT. SECOND. ^ORT. GIVEN 0-107. TOP KILL 

BP(IJ=f OR 10)=EXPECT. SFCOND. MORT. GIVEN li-257. TOP KILL 

BBtlt J=5 0<? 11.)=EXPECT. SECOND. MORT* GIVEN 26-502 TOP KILL 

l?R(ItJ=6 OR 12) -EXPECT. SECOND. MORT. GIVEN 51-997. TOP KILL 

BP ARRAY IS NOT OUTPUT BUT USED TO RECALCULATE THE 

AMORT OUTPUT MATRIX 



P^WIND 13 

RFWIND 1* 

RFAD(13 $ 7000M ( .AMORT ( I , J| , J- 1, 7) , 1=1 , 7 ) 

REAO(11,7001) f(RR(I f Jl t J= 1,6) , 1=1,2) 

RFADlf ,7000 ) ( (AMORT (I,J),J=l,n,I = *,l<i 

R r AD(l^,7001) ((BR(I,J) ,J=7 f 12) ,1-1,2) 

7000 FORMAT 

7001 FORMAT( 



7713 



90 



91 



COMPUTE THF TOTAL SECONDARY MORTALITY FROM 
TTS TWO COMPONENTS 



12 



00 770 J-1. 

PRC1, J)=Pf 1 ,J)*RB(2,J1-CBB(1, J) 

SFT UP POST-OUTBREAK EXPECTATIONS FROM INPUT DATA 
DO 100 T=llt 

TRANSFORM CUMULATIVE TOP-KILL EXPECTATIONS TO 
EXPECTATIONS BY TO^-KILL CLASS 

DH 90 J=l,3 

JJ=8-J 

JJ1=JJ-1 

AMORT (I, JJI=AMGRT(I, JJ)-AMORT( I,JJ1) 

REDUCF TOP-KILL EXPECTATIONS RY DIRECT MORTALITY 

DO 91 J=2,7 

AMORT(T,J) = AMORTtI,J)*( 1 . -A ^OR T { T, 1 ) ) 

IB=13 

IFCI.LT.8) IR^7 



115 - 

0116 - 

0117 - 

ana - 

0119 - 
10120 - 



DO 92 J 
JJ-10-J 



92 



c 

0122 -C 


c 



Q12<4 

10125 

10126 

)G12? 

)0128 

30120 

00130 

00131 

00132 

00133 

0013^ 

00135 

00136 

00137 

00138 

00139 

001M 

flOlM 

001^2 

001^3 



001^5 



001<+7 
001^8 
001^9 
0015C 
00151 
00152 
00153 
0015** 
00156 
00156 
00157 



- 93 



- 100 
-C 

-C 

-c 

- 777 

-C 
-C 

-c 



-c 
-c 

-c 

-c 

-c 
-c 
-c 
-c 
-c 
-c 



-c 
-c 

-c 
- 1 



00160 -C 

00161 -C 

00162 -C, 
00162 - '< 



r 

1,1 W J- *J* 

00166 -C 

00167 -C 
0016* -C 
0016? -C 

00170 -C 

00171 -C 

00172 -C 

00173 -C 

00175 -C 

00176 -C 

00177 -C 

00178 -C 

00179 -C 

00180 -C 

00181 -C 

00182 -C 
00182 -C 
0018^ -C 

00185 -C 

00186 -C 

001 87 -C 

00188 -C 
01318 -C 

00190 -C 

00191 -C 



CONTINUF 

REWIND INPUT LUNS 

00 777 I = 



PARAKTERS FOR GROWN DEFOLIATION EQUATION 



P3=-5981 



n 

.0 



INITIALIZE OUTPUT FORMAT CONTROL INPTCES 



NEW 

i = n 

IFLAG-O 



PflGFt RF.SET LINE COLNTER5 



BEGI 



N NF.W LINE* INCREMENT LINE COUNTERS 









__-.-.- 




7(6) IN SO 




ICC 1,1) 
iCCIf 2) 
IC(I,3> 


HOST SP^CInS 
Tf<FE CLASS NO. 
LINE INDICATOR FOR 


J, COUNTER 
AMORT TABLE 


FOR P 


IjjilS 


NOMINAL 7. NEW FOL. 
NOMINAL TOTAL FOL. 
ACTUAL NEW FOL. 


ZC21 IN SO 
7(75) IN SO 
X<5) IN SO 

X<6) IN SO 




C<T,5) 


'/.CROWN DEFOLIATION 






win) 


WEIGHTING FACTOR 


7( c ) IN SO 




A|I,1> 


NO. OF TREES 
NO. OF VIAELE EGGS 


7(if) IN SO 

xm IN so 


IN SO 


A (I I 3) 

Atije>) 

Atlf 7) 

ACT, 8) 

A t T Q \ 


TOTAL FOLIAGE 
ACTUAL 7. NEW FOL. 
7. DEFOLIATION 
EGGS LAIO 
REDIST. FGGS 
NO. OF VIAf?LF EGGS 
THTfii FOI IAGE 


F^OH Y(?) 
Y(3) IN SI 
G(2) IN SO 
XC7) IN SO 
X(5) *X(6) 


IN SI 
IN SO 



U U 1 ^ ~ 

00196 
00197 
00198 
00199 
00200 
OC201 
00202 
00203 
0020^ 
00205 
00206 
00207 
00208 
00209 
QC21C 
00211 
00212 
00212 
00214 
00215 
00216 
00217 
00218 
00219 
00220 
00221 
00222 
00222 



-C 
-C 
-C 
-C 
-C 
G 
-C 

C 

-C 

-C 
-C 
-0 



A(I,l*f) 

AtTtlSl 
A(I,16) 



Ad, 18) 



A(I,?11 
ACT, 22) 
A<I,?3) 
AT,24) 



NC. OF VIABLE EGGS 
TOTAL FOLIAGE 
ACTUAL X NEW FOL. 
7. DEFOLIATION 
EGGS LAID 
PEDIST, EGGS 
NC. OF VIABLE EGGS 
TOTAL FOLIAGF 
ACTUAL '< NEW FOL. 
7. DEFOLIATION 
EGGS LAID 
REOIST. EGGS 



X(7) IN SO 

X C5 ) +X (6) IN SO 

FROM YC2) IN SI 
YC3) IN SI 
G(2) IN SO 
XC7) IN SO 
X(6I+X(51 IN SO 



F?OM 

Y(3) 
G(2) 



A(I,26) V OF TREE NOT DEFOLIATED 

RFAD(l,inni) ICM,1) ,A( 1,1) f WT tl),0(l, 
C A (1, 2) , 1C (I, 2) 
1001 FORMAT (I2,7F11.5,I3) 



1,2) 



Y(2) IN 

IN SI 
IN SO 



SI 



00225 
00226 
00227 
00228 



RFAO<2,1002) C(I,M,A(I,6 

- 1002 FORMAT (12X,2F12.6) 

- READ(3t 1003) A (I ^7) 

- 1003 FORMAT(1F12.6) 

- lOQff FORMAT* (if6X,3Flll5 f $X) 

A(lJlO)" = C (I,?) MQQ./A (1,9) 
READ(5,1C02) C ( I ,*f) , A (I , 12) 
ACT, 11) =100. 0-C(I, ^1*100. 0/CCI, 21 
PFAD(6,10H3) All, 13) 



00230 
00231 
00232 
0023? 
0023*4 



00236 
00237 
00238 
00239 
002<*C 



QQ2M 
002^5 
002^6 
00?i*7 
002M 
002^9 
00250 
00251 
00252 
00253 
00254 



-C 
-C 

-C 



-C 
-C 

-C 

-C 
-C 
-C 
-C 

-C 



A ( I,15)=C (1,3) ^C (I,U) 

A(I,16)=C(I,3)*100./A(I,15) 

READ(8,1002) C(I,f),A(I^l) 

A (I, 17) =10fl.O-C(I,U)*100.0/C(I,2) 

PFAn(9tl003) A (T ,19) 

PFAD U 0,1 00^1 C(I,3),C(T,*),A(I,20) 

A(ll2?)=C(ll3)*10nIo/A(I,21) 
PFAD (11,1.002) C(I,f) ,A <I,2M 
A(I,23)=100.0-C(I,^)*10C.O/C(I,?) 
RFAD (12,1003-) A(I,?S) 

CALCULATE MAXIMUM 7. OF FOLIAGE WHTCH SUPVIVFC I 
X1=PMAX (A (1,5! ,A(I, 11) ) 



IT FI&ST LINE, GO WRITE PAGE HEADER 
IF (I.FQ.l) GO TO <* 

CHFCK FOR NEW HOST SPFCTFS OR DIFFERFNT NOMINAL 
FOLIAGE CCMPLPMENT. IF DIFFEPFNT, SET FLAG ANC 
START NEW fAGE 

IP ( IP(T-1, 1) .NF.ICd, D) J 



00255 - 

00256 - 

00257 - 
0025* - 

00259 -C 

00260 -C 
OC261 -C 

00262 - 

00263 -C 



-C 
- ^ 



- 15 



OC?65 
00?66 
00267 
00268 
00269 
0027C 
00271 



I NCI ~IF IX (10n.*C (I- 1,1) + 0.5) 
lNC2=IFTX(100.*C(I,l)+.5) 

IF (TNni.NE.INC2) JMPFLAG-1 
TF< JMpFLAG.EQ.l) GO TO 3 

IF LAST LIN ON PAGF, SET IFLAG=1? WRITE LINE. 
!*=( I.EQ.16) TFLAG = 1 



CRO^N nE p OLlATTON ANC FIND ASSOCTATEP L! 
OF MOPTALI T Y AND TCP-KILL TA^LE (A^CRT) 

CONTTN'.IE 

IF(A(I,17).GT.<=5.) GO TO 15 



. 0311 



GO TO 16 



00275 





00276 





00277 


~ 


00278 


* 


00279 





0028Q 





00281 


-C 


00282 


-C 


00283 


-C 


0028<< 


-C 


00285 


* 


00286 


-> 


00287 


. 


00288 





00289 


~ 


00290 


-C 


00291 


-C 


00292 





00292 





0029*4 


-C 


00295 


-C 


00296 


-C 


00297 


- 3 


00298 


-. 


00299 





00300 





00301 





00302 


- 12 


00303 


-C 


0030*4 


-C 


00305 


-c 


00306 


- 18 


00307 





00308 





00309 


- 


00310 





00111 


-. 


00312 





00313 


* 


0031*4 


* 


00315 





00316 


-in 


OQ317 





00319 


^ 


00320 





00321. 


- 


00322 


_ 


0032? 


- 




- 210 


00325 





00326 


~> 


00327 


* 


OH328 





00329 





00330 


- 


00331 


- 


00332 





00^3? 





0033*4 


- 220 


00335 





00336 





00337 





00338 





00339 


. 


007^0 


. 


Q J 03^3 


: ?3o 


00^5 


- 


003^4? 





003*9 





QQ35C 





00351 





00352 


- 2^0 



TF {TGT-*5.*ANO.T.LF.65*> ICC 
IF (T.GT.65..AND.T.LE-85.) TC( 
IF IT.GT. P5. .AND.T.LE.95.) 1C (If 3) =5 



IF 
IF 



(T^GT. 95. ANO.T.L:99.5) IC(I f 3) =6 
ICCIf3)=7 



IFCTCCIf D.EQ.2) 1C ( I , 3 ) = IC ( I, 3) +7 



WRITE LIN(I1 FOR TABLE 1 
ON UPPEP HALF 0^ PAGE 



AND ? 



WRITE (15 * 20 02) 1C (I f 2) 



(I f JJ) f J J= 



3) 



,3002) IC(If2)fACI f l),WT(I)fACIf2)fACI51f/J 
CACTf 17) ,A(If 23)^ CCTf 5),<AMORT(ICKt JJ) t JJ=lf*) 



CHFCK FOR FULL PAGE ANC EOO 

IF (IFLAG.EQ.1I GO TO 18 

IF (INDEX. LT.INOFXLTH1 GO TO 2 



COUPLE SPACE TO BOTTOM HALF OF PAGE 



CONTINUE 
IJ=16-I 

no 12 jj=ifij 

WRITE (15,2002) 
WPITF (16,3002) 
CONTINUF 

FINISH PRESENT PAGE HALF 

WRITE (15,2003) 

WRITE (16,3003) 

IX=I-JMPFLAG 

IF <IFLAG.EQ-1) TX=15 

DO 10 II=lfIX 

WRITE ^1^,2302) 1C ( II , 2 ) , A ( I I 



WRITE (16^300*) IC<IIf2),A(II,l) f WT(II) f A (II f 2) , 
HCAMORTdCK, JJ) , JJ=5f 9) , 
CA(II,2B) , AtII,26) ,Af IIf26) 

CONTINUE 

00 200 Kl=lf* 

W^ITE(lflflOOOO) J, (TSPECIESIJJ) t JJ=IS f IM) ,C(lf 1) 
C 



K3=K1 



DO 210 K?=1,IX 
WPITFC18, 10100) 
CONTINUE 
ENOFTLE IP 
CONTINUF 



TIME^DATE 
WRITF (22^21000) 
00 220 Kl^lflX 
ICK=IC(Klf 3) 
WRITE<22t ?1100) 

'CONTINUE 

2? 



C 



J, (ISPFCIES( JJ) f JJ=ISf IH) 



f 1) fCC 



IC(K1f2)fA(Klf2)fA(Klf5lfA(Klfll)f 
) , AMORT (ICKf 1) 



(TSPECIES(JJ) f JJ=IS, 



,C(1, 1) 



WPITE(?3 f 22000) 

DO ?30 K1=1,IX 

ICK=ICCK1 f 3) 

WPIT<23 f ?2100) IC(Klf2)fA(Kif2)fC(Klf5)f 

< AMOPT < ICK,KJJ) , KJJ-1 t <+ ) 
C N T I NU E 
ENDFILE~23 
WRTTF2^f?00-00) J, ( I SPEC IES ( JJ ) ,JJ-TSfIH) ,C (If 1) 



WPITE (2^,23900) 
00 ?/4C <!=!, IX 
ICK=IC(Klf 3) 

WPITr(?^,?3100) IClf2)fACKlf2) f AMORT (ICK2) f 
C ( AHO^T ( ICKf JJ) f JJ = 5, 9 ) 
CONTINUE 



re . 



- 20 
-C 
-C 

-c 

-C 

fj 

-c 

- 19 



- 7 



IF NO HORE DATA, CHECK JMPFLAG FOR RETURN 
IF MORE DATA,, CHFCK TFLAG FOR FULL PAGE. 



LT.TNDEXLTM) GO TO 20 
TFUMPFLAG.NE. 1) RETURN 
GO TO 19 
IF! IFLAG.EQ- 1) GO TO 1 



*** LINEII) 

REPLACE 
REPLACE 



CONTINUE 

DO 7 11=1,26 

A (I, II) =A (I II) 



1C {1,1) =IC!T f 1) 
10(1,2) =IP<I,2) 
C(lt l)=CCTt 1) 



IS FIRST LINE OF NEXT PAGE, 
A(1,J) WITH A(I t J) f J=l,26 
1C (1,K) AND C (ItK) WITH 
AND C(IfK)fK=l*2 



-C 

-c 
- I* 



WRITE OUT HEADER AND BEGIN NEW s>AGE 

rONTTNUE 

J=J + 1 

1 = 1 

JMpFLAG=0 

IS = 1 

IF (IC(1, 1) .E0.2) TS = <4 

I" = IS* 2 

CPMP=100 .-0(1,1) 

WRITF Cl c ,2001) Jf 

WRITE (16,3001) J, 

GO TO 5 



t CATE, (ISPECIFS( 



,C ( 



-C 

-C ***** TQRMAT STATEMENTS ***** 



2,*t INPUT-OUTPUT VARIABL 
C YC L* , ^2X , 2 ( 2X ,A8) 9 // 



- 2H01 FORMAT(#l#t9X, XT4BLE 
C* VALUES OVFR T HF - 
THE *,3A*f, 

*,Fe.. ~S* TREF 
/3QX,*TNPUTS*21Xi*OUT p UTS#j21Xt*TNPUTS**21X 

* *129(*-*)/#0 TREE NUMER OF *,?(* VIA 

TOTAL PFRCFNT CEFOLIA- VIABLE REDIS- *>, 

CLASS STEMS t^(t EGGS FOLIAGE NEW*,6 

Ff-GS TRI r UJTED t)/t NO . *, 1*X , 2 ( 1 OX ,*BIOM A 
C#FOLTACE*ii3X,*LATD EGGS *)/?tSX,2(* (NO.) t 

C*(GPAMS)~ (X) (X) (NO.) (NO.) *)/* *, 

C1 29(#*) ) 



C^ PHASE 

C^O 

r*i" 

C/i 



- 2003 

-C 
-C 



FORMAT*l* f 9X*TABLF 2*^,12,*? OUTBREAK DEFOLIATION 

THF FXPLCTFO MORTALITY AND GROWTH REDUCTION*, 

*FOP THF * * 

?,* TPEE CLASSES*///, 12 C 5(*-*)/ 
HFTGHT INITIAL : MODFL 

pc?r >POPTlON OF TREES RECEIV1 
FACTOP VIABLE * PERCENT OEFOLIATI 



^r TPF C NUMER 

C17X, iPE PC r NT OF 
Z* Ci ^S 
P^ BY PHA 
r* NO, 
C*OI?FCT 



(/-/ ) ) 



OF 

, flX f 

TP^ES " PFR 

SETCNOARY 
PhflSf I 5 
i HOPTALITY 



EGGS tt,3(#-#) t # TOTAL 
TOP ONLY GROWTH^/ 
I*,7X,*II/,7X,#TIT*,eX,*IV 
MORTALITY KILL REDUCTIC 



-'?002 



*, I3*3Xt ? 



,X,F7. 3, 3X f 

\y * 



^003 FORMAT(7tO#, 120 <#-#)// TREE 

^ex,/pROPOPTIo^ OF T^EES RY^^TX,*: 

C/DIAMFTER GROWTH IN*/ 

C* CLASS OF FACTOR VIABLE 



NUMBER WFIGHT INITIAL 
,*: GROWTH IN SURVIVING 



? PFPCENT REDUCTION 



C*CROWN HEIGHT S TP^S WITH NO TOP KILL t TPEES WITH TOP * 
C* NO. TREES PER EGGS f *,33<*-*) f * : *, 

C*PERCENT OF NOHINAL : *,19(*-*)/ 
1437 - ni8X?*TRE H/iSF T S 00 60 17*5 37,5 70,0 ^ 

343P - C*DTAMFTFP HEIGHT 5 PERCENT OF NONINAL*/X, 1.20 (*-*)) 

3439 - 30134- FORMAT<*0 * * 13, 3XF6 . 1 1 3X,F5 . 1 , 3X , F5<, 2 , 3X f <* (F* .*4<> X) f X 9 



1000(3 FORMAT it TABLE 1,#,T2,#1 INPUT-OUTPUT VARIABLE VALUES * 
C *OVFR THE OUTeREAK^/^X^^CYCLE FOR THE * , 3A**,F5 , 2, 
#/# f f'5.2 f * TREE CLASSES.</57X t A8/35X, -#PHASE # f A'3 f 
C 13X t A5/lX f 13(#-5t)t #f -------- INPUTS ------- 1#, 

C #[ ------- OUTPUTS -------- ]*/* TREE NUNRER INITIAL TC 

C # PERCENT DFOLIA VIABLE RECIS-#/'# CLASS OF VI/J 

C ^^LF FOLIAGE NEW -TION EGGS TRIGUTED */* NO, 

C * STEMS EGGS BIOMASS FOLIAGE ^n 1 IX^LAIO 
1450 - C 16Xt*(NO.) CGRAHS) (XI # ,6X f # (X) (MO*) 

)451 - C iX,64(^*#)) 

1452 -10100 FORHAT ( lX f 13 , 2X ,F7. 2,lX f 6( IX, F7.3) ) 

1452 -20000 FORMAT!* TABLE 2#fI2,jM OUTBRFAK DEFCLTATICM SUMMARY ^^ 
Jf54 - C#AND THE ASSOCIATFO#/5X, ^MORTALITY PROBABILITIES FOR THE 

1455 - C,3A4 f F5.2 f #/< f F5 ,2/5X ,*TREE CLASSES* #, 22X, 2 ( 2X ,A8) /) 

)456 -?100U FORMAK* ------ A: DEFOLIATION SUMMARY *,35 (#*#)// TREE < 

1457 - C* INITIALS, 5X,<RPANCH PEPCFNT PERCENT CF EXPECTATK 

!458 - C/# CLASS VIAPLE DFFOLIATION HY PHASED, 7X, *TRF#, 8Xt #OF# 

1459- C* NO. EGGS *,22(*-#),* TOTALLY CIPECT*/7X, 

1460 - C^PHASF II II III DFFOLIATFD MORTALITY*, 

1461 - C/* *,6<* <*-#* 1 

1462 -21100 FORMAT(2X,I3-t 1X,4(1X,F7,3) t 2X,F 7.3 5X, F7.*> 

1463 -??SCin FORMAT!* ------ MORTALITY PROBABILITIES *,31(*-#)/ f 

1^6^ - C* TREE INITIAL PEPCFNT*, 1 3X,*XPFCTATION OF*/* CLASS* 

I46S - C,* VIARLF OF TREE #,4i(*-*)/* NO. EGGS OEFOL-*, 

466 - C* DIRECT SECONDARY TOP NO TOP K ILL^/7X, *PHA SF I*, 

1467 - c* IATED MORTALITY MOMALITY KILL DAMAGE*// *,6*(*-* 



469 -23000 FORMAT C* ------ G? TOP KILL PROBABILITIES *^lt*-*1/, 

I47G - C* TRFE INITIAL*/* CLASS VI ABLE* 1 1 <X, *EXPECTA T ION OF TOP 

471 - C* KILL*/* NO. EGGS *,5Q (*-*) /7X, *PHASE I ANY L 

**7? - C*EADER 1-107. 11-2 I 5X 26-50X 5J-90X */* *,6M*-*)> 

473 -73100 FORMATC IX ,1 3 , 3X t F7, 3 , 2X , 6 <F7.4,1X)> 
47i - END 

F CN INPUT UNIT 



00001 - SUBROUTINE TABLE32(B) 

00002 - DIMENSION B(l) 

00003 -C 

0000^ -C WRITE THE DESTRUCTION/CONSUMPTION RATIOS AND 

00005 -C INDIVIDUAL PIOMASS GROWTH RATES ON TABLE 3, LUN 1 

00006 -C 

00007 - WRITE (17,1500) (B( J ) , J=7 , 12 ) 

QOHOe - 1500 FORMAT** # , 2< ( #-#1 /# DESTRUCTION/CONSUMPTION*/ 
QOOQC - H* #, 10X, RATIOS*/* #,2M#-*)/#n#,SX 9 #NEW FOLIAGES, 

00010 - C8X, #1,11, III QF,GF#,5X,6(F6.3,X11 

00011 - WRITE (17,1550) (B ( J ) , J=13, 1 8) 

00012 - 15^0 FORMAT(#0#, 8X,#Otn FOLI AGE#, 8X , #1, II , 1 1 1 tF,GF# f t 5X, 

00013 - r6CF6.3,X1 ) 

0001*4 - WPITE (17,16QQ1 8(22) ,?(221 ,8(221,3(231 ,<M241,8(2*1 
00015 - 1600 FORMAT(* t , 2^ ( #-#) / #0*, *5X, t IN DIVI OUAL ift BIOMAS's^ 

OOOIE - C* GROWTH RATF*,5X, #1,11, III CF,GF# ,5X f 6 (F6 ^, X) ) 

00017 - RFTURN 

0001P - PND 

00019 -C 

000 20 -C 

00021 -C 

00022 -C 
0002? -C 

OQn2<* - SUBROUTINE S2 ( X, KP,.F040<f) 

0002S - DIMENSION I BOUT ( 8) , IXOUT (8 .) , XOUT (4) 

0002 - DIMENSION X(l) 

00027 - EQUIVALENCE ( IBOUTf 1 ) ,X1) , ( I80UT( 3) , X2 ) , ( IEOUT (5) ,X3) 

OOH2^ - EQUIVALENCE f IEOUT7) f X41 , (IXOUTdl , XOUT( 111 , ( IXOUTC3) , 

00029 - EQUIVALENCE f IXOUT (5) ,XOUT ( 31 ) , (IXOUTf 71 ,XOUT) ) 

00030 -C 

00031 -C AT KP=0, WRITE CURRENT STATE VARIABLES ON LUNS 21 

00032 -C AT ALL OTHER TIME STEPS WRITE UPDATED 

00033 -C VALUES ONLY ON 21 AND 25 

00035 - Xl=X(1l*(FCl,l)-Ffl, 211*0. 10 

00036 - X? = X(2)+IF(2,2H-F (1,211*0. 10 

00037 - X3=X(31*F(3,3) *0.1Q 
0003? - X**=:X ( **) *F O^f *t* 0* 10 

00039 - IFCKP.EQ.01 WR ITE ( 21 , 100 0) X (11 ,X (21 ,X (31 , X (^1 

000^0 - W^ITE (21,1000) Xl,X2,X3,X<f 

000^1 - 1000 FORMAT CM ?X ,FH.71 ) 

000^2 - IF (KP.FQ.9) WRITE (21,1000) 

000^3- 00 1 T=l, <f 

000^ - 1 XOUT (I) =X (I) 

000t*5 - IFCKP.^Q.O) BUFFER OUT(25,11 ( IXCUT ( 11 , IX OUTI 81 1 

000^6 - CALL UNITST (25,IR?5) 

000^7 - IFCKP.NC.9) BUFFER OUT(25,11 ( I ROUT ( 11 , IPCU T(8 1 1 



000^9 - END 

EOF CN INPUT UNIT 



coutine l 



riowcnart 




Note: The numbers inside symbol FH are statement labels in the code for 
subroutine TABLE12. ^ 



Subroutine TABLE12 flowchart (continued) 



LINE I BEGINS NEW 
PAGE, TRANSFER 
TO LINE 1 STORAGE 
VARIABLES 
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APPENDIX C VERSION 3.1 VERIFICATION 
INPUT AND OUTPUT 

Page 
Inputs : 

LUN1 79 

LUN13 79 
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The value 0.333 for this parameter was used for validation simulation 
i updated value 0.030 should be substituted for any further use of this 
ameter . 



D 
1- 



X 
0. 



O 

ILU 



JDCDDZ 
UJMU*- 

arar 



coot- 
<ci 

MU 

> 

t 

*c 

M2 
-JO 



UJ 

a 



ui 

a iL 



o 

AMJ. 

- ._ .DO 
OCD03C 
UJMl J*-- 



UJ 
-ICOC 

WLU.Jw 



UJ(1- 
JIOtT( 

<t<t 



H 


03 


ro 




4 < 


3 


fs. 


N. 


3 CP 


CO 




o 


CM 


ro 


in c 


3 4- 


O3 




3 O 


00 




1-4 


CO 


4- f 


CM 


in 


ao C 


* in 






03 




N- 


5 <H 


ro 


K L 










CM 


CM 


j ro 


ro 


ro , 


r 4- 


H 


03 , 


-t ro 


4 


4" C 


3 


fs. 


N- fl 


D CP 


CO 


O i 


4 CO 


ro 


in c 


3 4- 


03 


CP t; 


3 


CO 


O C 


9 -H 


CM 


4- r 


- CM 


in 


eo C 


^ in 






OD 


4- 


N. 


5 T-J 


ro 


K t 


1 4- 








CM 


CM C 


j ro 


ro 


ro 




CO 


O C 


a CP 


CO 


CP * 


H 


j. 


ro < 


3 K 


CP 


CP 


o in 


ro 


oo * 


h o 


cr 


in 


^ CD, 


CO 




P cp 


cr 


CM C 


J 03 





C 


CM 




4 i 


N 


O3 


CM C 


* in 


CM 


in c 


3 CO 






H 


CM 


ro r 


"> 4- 


in 


03 C 


* CP 


K 


N- 1 


r\ CP 


ro 


cr v 


3 CO 


cr 


H * 


t 03 


CP 


CP t 


!> K 


03 


in 


r 4- 


ro 


CM C 


* N, 


CP 


CP < 


P CP 


CP 


cr c 


r cr 


<p 


CP 


3 oo 


CM 


CM 


CO 


CM 


CM C 


J CM 


CM 


CM C 


J CM 






































CM 


H 


HI CM 


4 


4" I 


"> 03 


03 


K e 


s co 


00 


H t 


O 4* 




h- < 


D CP 





CM f 


- -H 


CP 


CP I 


D OD 


4- 


r> c 


J -H 






^ CM 


CP 
CP 


CP ( 


P CP 


CP 


<p c 


* CP 


en 


eo 
CP C 


; o* 

> CP 


* 


CP C 


P CO 


03 


03 


r in 


4- 




, CO 




1-1 1 


O |s. 


H 




rv K 


CP 


ro - 


HI CO' 


^ 


(S. . 


t eo 


ro 


o r 


" **" 




O3 I 


n ro] 






<H 


CM 


CO C 


o ro 


ro 


4- || 




















! 


















i 


00 


CP < 


D IT> 


ro 


CO - 


H cr 


00 


03 ^ 


H Is-E 


eo 


ro I 


- in 


ro 


f 1 - ' 


H .4- 


CO 


OJ C 


is. 


4" 


4" < 


D fs. 


03 


I 


r\ cr 


ro 


CM C 


D ^ 


H 


N- , 


t CP 


Jr 


CM ( 


P 03 




CP C 


n o: 






H CM 


4- 


in i 




r- 


oa , 


H 4" 


CO 

oo 


CP 


D in 
- in 


ro 
ro 


K I 


is 


co 


^ 

03 , 
03 C 


H ^i 

^ rs. 

J IS, 


4- 


4" 


8 N. 


03 


O L 


"i cr 


ro 


CM C 


3 rs. 


^ 


N- . 


t CP 


4 


CM c 


n 03 


4 


CP q 


n eoi 




1 
1 


^ CM 


** 


in i 


03 


N. 


00 i 
- 1 


H 4"! 
1 ,rtj 


CP 
-H 


03 I 
CP < 

in ' 


o ro 
n oo 


in 

fs- 

in 


H r 

i 


- ro 


in 
cr 


in r 

<H U 


r is, 

? CP 1 






-t CM 


ro 


4 .J 


t in 


in 


N C 




O 


t 


9 


o 


C 


D 


CD 


O ti 


3 CD: 


O 





a o 


O i 


o 5 


> o 


, 


CD Cj 


3 0- 


a 


o c 


3 


o 




3 


CD 


fi 


3 CD; 





C 


3 CD 





O ( 


3 O 


C3 


tt 


3 0- 


ro 


ro r 


* ^ 


ro 


ro rf 


> ro 


PO 


ro r* 


i ro; 








3 


o 


o er 


> o 


O : 


o et 


3 O| 




O tt 


3 


o 


o <i 




CD ; 


O e 


J 0| 





CD <fc 


3 


o 


$ 


3 


O ; 


CD C 




o 


o <e 


J O 


CD 


O cp 


, 


o i 


efc 


> CZDj 


CD 


tfc 


) O 


o 


CD < 


J O 


CD 


o c 


J C3 


o 


o i 


> CD, 


O ' 


o < 


> o 






) a'< 


O 


o i 


7 O 


O ' 


cj; 


J, ^5). 


O 


o cc 


J 0- 


v-4 


in c 




o , 


UN C, 


j in. 


o 


o ct 


> CD 1 . 




t- 


^ CM 


ro 


ro d 


r 4 1 - 


in 


VJU OC 


> CD; 






i 












^ 


O 


Ci <t 







o cp 


> CD' 


LD 


C 


) 


O , 

o 


O (t 

<t 


> 0' 

> o! 





o cp 

o or 


> CD 


O ' 


o c 

o c 


> CD 



or cr o; 

H -r< CM! 



CM CM CO CM 





L 


N 4" 


50 


r-. , 


4 4 1 


vO 


CP 


f tO 


H 




( 


r> K 




CM ( 






^ . 


" CV1 


O 




C 


9 J. 


CP 


n c 


3 O 









ca 




t 


^ 4" 


co 


rs. , 


4 Jf 


03 


CP 


h o 


rH 




c 


i% K 




CO ( 


r> <p 


H 






a 




c 


D 4 1 


CP 


ro c 


3 O 





cp ^ 


r M 







c 


3 CM 


in 


i 

N. i 


> 03 


H 


03 - 




CO 




c 


3 O 


a 


4- r 


5 ro 


O3 


in a 


J f<> 


eo 




\. 


3 





o ., 


H ro 


03 


H J 


K ^ 


CM 






ro 


03 


CO V 


3 4- 


JI 


H C 


ro 


in 


t- 








4- U 


> 03 


N. 


c 


> O3 




u 




















c/ 




CP 

4 CM 


4- 


4- U 


> CP 


4- 


H e 


9 W 


CM 


I 


r 


- in 


CP 


CP t\ 


J in 


O 


N- C 


3 CP 


co 




c 


^ eo 


03 


03 


H is- 


03 


in c 


3 K 


in 




c 


J CO 


CO 


ro . 


f 4 


in 


N- C 


1 4" 


09 




a 


3 CM 


ro 


03 C 


* in 


CM 


4- 


-H 


^ 




e 


p 4 


CO 


CP r 






cp c 


> CJ 


1-1 




L 


N CP 


oo 


in < 


3 CO 


ro 


N ^ 


r ro 


PS, 




< 


N. 


in 


4* i 


3 N" 


o 


O 


< Is. 


4" 




C 


> CP 


CP 


in . 


t ro 


ro 


CM t 


' w 


CM 




C 


*> 4- 


CO 


ro < 


r. 03 


ro 


CM r^ 


ro 


03 




C 


J 4* 


00 


C 


D 


N. 


03 r 


- N 


O 




t 


^ vO 


CM 


4" c 


3 -H 


03 


4 j 


C CM 


o 






CM 


in 


N. v 


D 03 


in 


in c 


J -4 






I 




ro 


4- C 


K 




CM 


in 






! 


CJ 
"\ 03 


in 

CM 


in c 
to i 


n o 


CO 


4- U 


4 CO 


rr> 

CD 




r 


^ K, 


in 


CP C 


3 O 


rs. 


03 ^ 


CO 








H 


ro 


4- . 


I- 4 1 


ro 


ro 








i 


> O3 

CM 


ro 
in 


4- C 

in c 


o N 
r> o 


H 


CM ^ 


in 

CO 


N- 

ro 




L 


"N vX) 


CM 


ro L 


*v N 


00 


JT L 


4 


CD 




r 


3 N 


in 


cr- c 


3 O 


fs> 


03 v 


CO 


* 






H 


ro 


4- , 


r 4- 


ro 


ro i 








i 


> ro 


03 




r> in 


CM 


CM C 





ro 




-4 


r- co 




O3 C 


3 ro 


4- 


03 C 


<p 


ur\ 


H, 


c 


j ro 


N 


IT* n 


3 


03 


4* C 







E 




rH 


CO 


03 r 


- OO 


03 


CP e 


eo 


CP 


uj 




















a 


v 

I 


3 in 

' St' 


CD 


CM C 

S S 


J CP 
* 00 
3 O 


ro 


-H r 

03 U 

CC U 


00 
K 
03 


cp 

O3 

CO 




a 


1*1 cp' 


CP 




! N 


CO 


cr i 


03 


CJ 


i 
1 


' 


J CMj 


CM 


CO P 


J CM 


CM 


CM r 


ro 


vO 




u 


^ N. 


ro 


OD ej 


> 00 


O3 


cr c 


cr 


CO 




VJ 


3 CM 


in 


U 


> vD 


03 


in c 


oo 


N. 




Of 


> ro 


03 


ro fl 


> flO 


CO 


in c 


H 








<! 


i cr 


CP 

i 


CP 3 


> CO 

i cr 


CO 


K f 


in 
in 


O3 
CM 




j 
it 


J CM 

3 03, 


i 

CM 


5 r! 


r CM 


in 

fs. 


vO J 
CP J 


. cr 


ro 

CD 




r\ 


J TH 


ro , 


O!> C 


s rv 


o 


O <vj 


4 CO 


03 




rr 


) OJ 
-H 


CM : 

ro 


03 C 


3 cr 

cj 


H 


u> J 
CM H 


> in 


ro 

CO 
























c 
c 


> cr> 


CD 
CD 


C 

CD C 


J 
3 


O 
O 


o q 
o q 


; z 









c 


3 CD 


CD 


O C 


1 


o 


o q 


9 



























rf 


4" 


in ; 


o: n 


CO 


cr ! 


O i- 


* CM; 


ro 






-H 


vHI ' 


<H T 


i -H| 


H 


CM Ci 


j CM) 


CM 


















! 














i 




i 











t 




1 


















































g 






















































1 


X 


I 


















































f 




' t 


















































t 


3t 


> 9 


9 C 





o 


3 


o 





3 





1-4 
































1 


O 






o 







o 










O 






i 


























1 


a: 




9 C 


o 


o 




o 


o 


9 


C3 


ro 






i 


























\ 




r> a 


cr 


cr> 


n 


en 


cr 


r 




































1 


> 


























i 
l 


























1 


_j 






















































z 




























j i 
























I 





























I Z 


Jt u 
























1 




























{ >H 


4t d 
























1 




























1 


f 1 Z 
























i a. 


























1 X 


H- 
























I O 

i (- 




r\ ir 


in 


in 


> i 


in 


in 


> i 


in 


CVJ 






1 3C 

J O 


3 1 O 


a 


cr 


00 


OC 


CP 


a 


J 


w 


N. 






i 

i 




























1 


: f U- 
- o 


I a 





a 


O v 


ro 


in 




r*; 


UD 






> 




























i or 


i I 


i a 


cr 


CO 


n -; 


in 


in 


n i 


in 








i 




























1 UJ 


C 1 - 
























i 































1 Z 
























i < 





























1 UJC 


"> 1 U 
























i 


tfl- 




* ^ 


^ 


^ 


<. ^ 


^ 


_jf 


^ 


^. 


CJ 






t XI 


J 1 O 
























i < 


3>, 




N I 


in 


in 


N If 




in 


^ l 


in 


in 






t <CL 


1 Ctf 
























i: 


Z 








jf- 




^ 
















^ 1 UJ 
























| 


D- 

sn 




1 C 


o 


o 


3 e 


o 


o 


3 





o 






1 01 


. t a 
























1 UJC 


























1 


MM 
























! ( 


/>S 


























1 
1 . 
























1 


, 




























1 . 

Ot- 


J 






















1 


i 


1 


























25 


<f,<C 






















I 


i 


H 






o 




3 a 


o 





3 a 





a 






1 t-i 
























1 


O-J 

) U-la 


























I > 


.MH- 


1 oc 


cr 


00 


JO QC 


cr 


o 


J 


ro 


is. 




1 


l ! 


Xt- 
HO 


























^>> 
1 Of 


-00 


^ 


o 





4 


ro 


^ 


r> cr 


vO 


uo 


Ul 


1 


> OC 


























1 OC 

l to; 


UJ 

U-X 


c 


0> 


00 


n .j 


in 


in 


n u 


in 


ro 


































i 























(X) 




M 


* 


























i z: 






















-J 


u. 




C 


























i * 


za; 




















o 

UJ 

UJ 


c 
- 


vu 
,J- 

AJJ< 




C 


o 

a 
o 


ro 

' 


i : 


cr 




< 


D CJ 
n N. 
3 t3 





a 

ro 






i x: 

1 2( 


LULU 
01- 


i c 


ro 


00 


OC 

> vC 


s 




in 


; ^ 


ro 

ro 


K 


Of 

H- 


c 


fj-, 

-OC 


















vO 


OO 
CJ 






1 Ol 
1 OA 
I OC 


JJX 

a<c 


c 


o 

cr 


a 
<r 


T* v 


ro 
in 


in 


n u: 


in 


ro 




a 


>. 


-u 


























I t 


o 




















ep 
o 


K 




! 




























1 * < 






















o 

a 











N. 

ro 


jr 

N. 


eo i 

OC ' 


VJ C 


00 
IS. 


CM 


CP 
CJ 


H 
O 









| f. 


1 O 





00 


00 


BO C 


oo 


00 


3O OC 


00 


00 


in 










rr 








s. 


IS, 


ro 





^ 






f : 


t 


c 





o 


c 


o 


o 


3 C 


o 


o 


































1 C 


1 


c 


o 


o 




o 




> C 


o 


Q 








5 










r 1 ci 




cvj r 




a 


cr 






l > 


4 1 I s - 










* 










v 




UJ 1 


K 








H ' 


VJ CV 


-< 


H 


H 




,^- 






1 L 


J 1 
















* 


* 






(/) 






























: i 




















o 
































1 >-" 


i 




















o 






X 1 

a i 

i 




c 

VJO 


<T 
CJ 


% i 


c: 
ur 


o 

LT> 


ro 
JT 


ro 

UJ 


00 
CVJ 
OD 


in 

H 






i co; 
i : 

1 UJC 


> 1 lf^ 
> 1 


r. 

c 
c 


oo 

ro 





00 
U 


SO QC 

=> c; 

=> G 


CO 

ro 




ro 

o 
o 


C 

c 


OO 

ro 
o 

o 


OO 

ro 



o 


ro 






- 


























1 UJC 


> i ro 
























i 


t- 




<T 


cr 


O 






O 1 




CO 


H 
































> 


k- 










cr 


4- 




J- 


vO 








!* 




















< 




-4 




1 






















\*s 


4 1 

i in 


cv 


N- 

co 


CJ 


VJ CV 


CVJ 


N. 

CM 


N 


N. 

CM 


M 


U- 


j 


H 




| '; 


00 

cv 





ro < 


00 

CT 


uo 

CJ 


in 


s 


CD 

jr 


i 
H 

oo 






i o; 

1 C 

1 ZK 


l 

> i N- 

4 1 1-4 


c 

c 









t: 
c: 


UJ 



CD 



C 

c 




o 


O 




o 
z 


< 


U 




[ ' 


*" 


(VJ 


in ) 


00 


o 


H < 


^ 


in 


C\f 




1 


















* 


<t 
cxr 


a 

00 


J 


1- 




N" 


vD 


CJ i 
H 


^ 


in 

CO 


CO 

cj r 


00 

ro 


r*. 
in 


oo 




i -: 
i tvc 

) OLL 

i 0. ck 


1 

1 
1 * 


rr 
cv 


ro 

CVJJ 


ro 

CM 


f. 
cv 


ro 

CJ 


ro 

CM 


cv 


ro 

CM 


ro 

CM 


UJ 


u 








1 


^ 


oC 


in f 


"> tv 


1 

O i ^t < 


vO 


H 


vD 






{gli" 


^ 












c 


o 





X 

~ 


o 

X 


1 




i ; 


Tj 


CM 

00 


in * 


3 (T 


H 


CJ f 





CJ 
JD 


CJ 






i t 


I 
| 






















u. 






f- 








i-l IV 


J CV 


ro 


ro 


* 




or 






i c: 

I C-, 

< i 

1 C. 


I o 
1 

| 

1 


c 


cr 

ro 
o 


ro 




f. 


ro 
o 


-3- 

ro 
o 


<* 


ro 
o 


*7 
O 




M 





, 
























; j ^J 




























a 


O 


O IS 


> O 


o 


a < 


a 


a 


o 




































< 






O i = 


j c: 


o 


o 


o 


o 






























< 

M 



M 


CO 
OU 
Of 




in 




o 


o <: 

CJ l^ 1 


) IT 


o 


IT 
-* 1 





o 

OP 


a 






! Ser 


OLU 


3 

H in 

1 - 


o 


5 



UPv ; 


C3 

c: 

t 


o 
o 


o 




Z 




3 

Q 




j 


i 










i 


i!^- 


1 


T-l 


H ro ; j- it ijr 


Vf 


OC' 




X 


CVUI 
UJUI 
0-QT 


<: 

f 


! ! 

> C3J' O O $3 Q 

W i* - 

4 -ri i-< | -^ 4-4 ^ 


o i o js a! 
; *| 

YH -H |H ^i> 


o 


H 


i i ; i 

1 MOCLfi'il 


CJ 0,0 k. 

.p- 


' 1 

4 1 4 ' 


o 




; 




T i : 


' ' , 1 LU- H- : ,' : 1 
























; ! ; ! : i i' : I ' i 




a: 

UJ 
CCl 

: 


CO 
UJ 
UJ 
Cf 


i 


f 

f 


f . ' i : 1 : || ' . ; j 


1 : i ' i 






D 





1 i ! j p | j . . j ^jTJl jj i i i ' ; 








; i i 


, 


' :i z ii ; = i 


1 UJ 
i Ul 


o 


+, cvj m 


* T a, K, ^r J ^! N ! :L,S. lS^-K,Luia,:KLJ,: rt 


1 z 

HO 


1 : i ; ' i * i ^ l ; i5* i 1 ! : ! ; i "V 




i ' ; : : ' .' : i : ! 



3 

z 

M 


TOP ONLY GROWTH 
KTi i PFnilCTTfJ M 


< 
1 
t 
i 

< 

L 


3 


1 

i 

N. |K f|0 CO 

in in PK o 
K ps. Jar ro 


eo 


"H V. 

CO C 

o r 

CO C 


D in 

9 CO 


H 

eo 


CO 


03 



O 


D 

r\ in 

-4 -H 


1 

i j 

i i 

CC> CO (W O 

in in ^p co 

v4 <H 4j CO 


CO 
eo 


in L 

o c 

C 


TV o 

P CO 
CO 


in 



o 


in 

CP 
CO 


IF WEIGHT INITIAL I MODEL 9RANCHf PERCENT OF 8 PROPORTION OF TREES RECE 


IS PER EGGS t TOTALLY f DIRECT SECONDARY 
Tf?rfr PHASF T ! T TT TTT TW nFFniTATFfl f MflOTAaTTY MnPTAlTTY 




VI CO 


co cj J3> in 


CP 


J- 


D in 




eo 

UD 




( 


>0 CO 


CO CJ 


CP 


<H 4* 


H 


( 

t 
( 
< 


3 

3 

9 
3 


o o b o 

<Jp CO 
<P CO 

i* 
cvj ro ijn cr 

J* vD 

CO P -rH 


ro 

N. 

vO 

03 


O 

N W co 
fN. CiO CM 
^ C|P o 
! o 

| 

CP tb K 

in <p > 


O 

" 

ro 
in 


o 

ro 

CO 

CP 

a 




< 


p co 


CP in co 

IJO N. 

in ! N. ijn u3 


O 
CP 


CP ip N. 


CP 

CO 
CO 
CO 


CO 

o 



o 

oc 


i 

I 

< 

t 

: 

< 




; 

d 


D 


o o <f"i ro 


03 


H 


p ^ 


r* ro 

t co 

S JS" 


' jf in 03 
03 co cr in 

J- 03 O fO 

j- i in *p -r-i 


H 

CM 

J" 
UJ 


H <JP ro 

CM O O 
VJ C3 U'. 


IT. 

ro 
in 


vO 

QO 

d 

o 
o 


vj ro 

H 

: o 
P CP 

r j- 


CO ! 0? fS. CO 
O CP <O CP 

or in m on 


CO 

N. 
CM 


CP cp ' 

H 

^ JT ro 
o cr in 


CP 

00 
CM 


CD 

H 

N. 
O 
(V 


P O3 

-4 cr 
-t in 


co ro in <H 
CP ; co h- r-~ 
- ro in -H 


CP 

ro 

5 


in co m 
j- in u3 

ro cr o 

03 in in 

ro dp -H 


CP 

ro 

in 


OC 
<T 

N. 

cr 


D 
J 


P 


H : CM PO *? Jf 

O O O Ol O 

o o o in; o 


in in N. 

C3 O CD 
O O LJ 
U" <P C3 


cr 

a 
c:- 

o 


-r-i 


O 

a 


O O O O' O 


d- in ur- 
ea o o 


o 


^ 












z 

LU0 
Q H 2 ' 

"p 


' 










T-4 --!- * -H T-' -H CM CO CM 


CM 


Cu 



. 
















1 
















1 \xl jj 

; i x M ! 
















IKCu X 
















3SO i 
















Ofr- 2: u3 00 


f\ $r 


CP 


U3 


uO ^ 





*"< 


n a 


! of co co 


n ^j 


<T 


CO 


ro c 




* 




1 1- N. j 03 


co d 


03 


CP 


*> if 




cn 


o 


i 1 Ct'l-t CP j OO 


- r< 


CO 


^ 


H 




H 




j t^* 2 
















i UJ(/J UJ 
















i ZLU O ! 
















! tbJ Cf | 
















| ! Si- a 
















: , co M ; ,oo 
















1 ! 
















: J ! 
















z vGc 
















! , >OLpl- 1 V0 i OO 


Lfv f* 


CP 


03 


X) ^} 


o 


T-i 


n c 


! : hH Ot X CO | CM 


m j 




CO 


*J C 








: Q, &M K i y0 


CM j 


03 


CP 


O IT 




CP 


P 


i : ZJO; LU CP ! on 
U02TL.X 


is. K 


CO 


H 


H 








i p 
















zx; 
















. wt-^. 
















j 0>- CO CO 


s s 


CP 

CP 


03 
of 


s 





3 


n c 


OLUJjjSC K ' O3 


ro J 


03 


CP 






cn 


D 


Ctf UtJO- I CP CO 


N. i*" 


CO 




H 




T-4 




: H' o 
































i 0,0 


G 


O 


O 


o a 


ro 





r r 


J o 








ec 


T-) 




T- 


xi 








H C 


O 






Oj Cj 















3 C 


, -! K 








i 








Uj' 
















5 x! 
















>. I . 
















onzr: o o 

3. in 


C 


S 


co 

ro 


n a 





CO 

ro 


3 C 


too; 




CO 


CO 






CO 




UJUC, i"^ 




o 


o 


s 




CD> 




LOO ro 










a 









KZ' 
















u. in 


J 






-, 








3 


O. N 
















Zt-n; -H 
















01- 
















^^' J- J- 


-* J 


in 


U3 


^ -4 


N. 


03 


3" N 


CfO C? O3 uD 




1 CM 










C LC 


OLU .* j- 


^ J 


( vH 


03 


in f 





U3 


*) a 




*a c 


> CM 


H 


f> c 





H 


D C 


a z' 
















LUT J- J- 

o o cr cr- 


CP 


: S 


UD 

N. 


JD (TV. 


U3 


03 


\JI 01 
CV 


ct or 03 


vD ^ 


i oo 


f"O 


"O if 


CO 


fO 


n c\ 


til CJ O O 


/=> C 


> T- 


CM 


H T- 


CD 


CO 


H C 


f ~ .. 
















-JUJ' -( 


C 


i O 


S 


C 


1 o 





P S 


5 ^^'5 


^H f\ 




M 


-* J 


in 


03 


pn o 

. 


> 1 o_ 
















1 O/ CD O 


'3 C 


> 





C 


> CD 


O 


b> c* 


HOQ. (V 
















Lu<l ~ 
















Ce 00 CJ 







- 




i O 


* 


o c, 


(/) 



































( 


1 
















































i 


UJ f 


























C 
























M 


























C 


ooc 




















( 

1 


M 
Irt 



























cv 


oini 

CM-H^i 




















i 


h~ 

M 






















OC 


OO 


























-J 






















oc 

in. 


Lno 
o en 






















SJ 




<* 


















































O 














































H 




X 














































M 























o 
































IH 

_J 


















in 


sci 






























<y 


















































o 

X 














































9 




vT> 






CM 

O 


O-H 

o 




BDin 
tMro 

00 


O- 
0< 


30 


OO 






















o 

o 




CD 


in 
cv 

v) 

a 




in 




j 








* 




* * 


* 


* * 






















CO 


ro 







CO 



















































in 




D 
O 




in 




o 

CM 




O-H 

o 




in in 
*>ro 

Oc3 


Oi 
O( 


pro 
30 


00 




















N. 
CM 


CO 


in 
<\ 
a: 




u 




"Bt 






































CM 


ro 






M 




D 








o 


oo 


ro^r 


oo 


3-r-l 


00 




















O 





u 












j. 




CM 

o 




oo 




4 CM 
30 


o 

H 




















K 


OD 


cO 
C3 




co 
ro 




J 
3 




: 


































CO 


w 






in 




O 




z 




O 
(\. 


oo 


vO-H 


00 


3O 


00 

OO 




















o 




in 


5 




^ 




3 




ro | 









00 




30 




* * 




















CM 


CO 

* 


. 




it 

M 




. 








* 




* * 


























03 


v> 






M 




^ 














































M 




< 
H 
U 
3 




CM i 




o 

CM 
O 


00 


roeo 

OCM 

oo 


oo 


30 


00 

r\o 

S.OC 




















O 

in 

CM 




in 

CM 


2 




10 

ro 










































vD 


uD 






CO 




- 









CM 


oo 


OCM 


OO 


3 CD 


00 

00 


























o 


5 




CO 
UD 
CM 




J 




H 




O 




00 




30 


00 




















- 


dr 


r- 




II 


O 














* 




* * 





















































<H--l 




















in 


in 






M 




g 

X 






















































/I 
























U-U. 




U. U_l 






U- 


u. 


u. 




O 




jJ 


j 


J 







oo 

9- 


oo 


00 


DO 


LL.LL. 




o 


Soc 


o 


>oo 


c 


OO( 


5 









> 


c. 




in 


< 


3D 








U, 

o 


oc 


DO 


00 


uU. 

DO 


OO 




u_ 
a 


JLU,U 

DOC 


u. 




.U.U. 
3OO 


u 

c 


U.U-1 

oot 


i 




u. 




U, 
O 


u. 

c 




i 


UJ 
Of 
















































M 













> 









































H- 










M 










MM 






























D 


- . 


f 
















r 






























O 


< 








t-t 










MM 










O 










M 


t-t 


l> 




** 




Q/ 




J{ 




M 




M 




H 


HM 






t I 


M 


^ 




M 






M 


hH 


>- 




U. 


X 


3 


', 


I 







MH- 


| H> 


MM 


H> 


* > 




M 


HM2 




t 


t~ 


MM 


: 




f. 









t-4 


H 


X. 










H 


MM 


M 


IM 








HMH 


1 






MM 


"H 




M 


14 


h- 
















H 










1|_ | 










M ( 










M 


f 1 


H- 




M 


a 


to 


. 






r 










> 








M 


*4M> 













*. 






(/ 


o 


* 


1 






M 










MM 








>4-IMM 










M 


I-H 


l_ 




Lu 


u. 


__ _ 














Q' 































h- 


jj 














O 












1 
















2 


x 














t 






UJ 








2" 










O 


UJ 


<c 






UJ 






t i 


















o 


M 




) | 


I 











<x 








N 


-i ; i 




u. 


O 






<c 




1 U.' 






ttT 






H ! 








M 


Q 




j 


Of 


I 










i '01 


X 


<u 


2 


U. 






O 

M 




O 

2! 


! a 

, N. 






j: 


i uj : i (/i ' 


T uj ; uj i nd 

-5 00-" 


I 


U. 
U, 


t* 




_ 


MO: UJ 




M 




1 (/) i 


/) i <I <I 


3 




3 


O 






IL. 


M 


O ! UJ ' t- 












zcn M ! M 








O 






11 


-4 


UL ! V) ! M 




00 


L 


LJU' 'Mi 


^ 


ovj -j ; _i i -jj 


c 


U 




UJ 
-J 

er 




(t 
O 

IU 


i 


0; < 
O </> 

<t M 




a' 
<t 




C/) 
UJ 

y 


O 

r 


<: i 3c i 

D| 'Of 

=>! ! UJ I 
! > i 


| u 

ii! 


C_*-i , O i O 
*vH- U. U. 

Da z ; o 


<tc 

OC 

hH 


(. 


M 
U. 
M 
O 


2 
O 

M 

h- 






13 




er o 




a. 






' O i : 


M ' UJ ' -I 


>(/ 






D 












1 


2 




: ! 


^ 


'20 






a 


CO 








y. 




! I 


O 




\ 


<< ' 


U i 


o< 


t 


(/ 


M 








-J 




1 i M 


i i i 2j i 


D I 


23 






&. 


1 


5 




1 i 


1 * 


: : 1 2 


21 I 


t- 




fc 


t) 


i j 


01 | i 


a 






a 






M 


! i 


1 i 





^ : 


jj 






; 


O 


i 1 


8 i i i ; o 




D i 




( 


UJ 


| 




; ! ' ' ! : ! ' ! 


; i 




t* 



/I O 

:z "2. 



Q 
>- 

'K U. 

/o c 



S. M 

8- U. 

u. 10 

2: ;E o 



M 
M 

M 

HH 
1 1 
















OOCtJ-HH 
o |e a 

tn ^ff-m oo 


o <D <i o o o e - U 

M 

-J 

: 


oSpH^ 


0-fl-Hi 







^ 








fojfcbinro 
rovocvjmeo 


H X 

! O 

j-^-trsf^uDtsjin ji af M 










2SPSSS5SSSS 


ruin 


eOHC 
OOi 


DUD.3- 

^^oln 


ississib 


SIXTTf TNSTfiR 
_^4W- -5-.CLIAGE- -- FOR OOUGL4S F 
III IV I 


rH CM fO j- in in *> N. o a 
CO -H CM Jh vD N. CO CT< C3 to K -H 

H <T> ioo vD ^f ro CM <H *H cy in CM 
o Q H ro in vO N. co cr o .* or 


inj-c 

OOCM, 
HCTfl 


3C\JI/% 
OC\JCT> 


oo oo <^ in o ip CM 


mcfior>fv.N 
OOCMkf OOCV. 


>rocM< 
irocot; 


oin in CM o 
nee dro 


CMOC 


raON. 


S.N.SI 


: 

^^! 


>H CM CM CM -H., 

1 
CMOJCT OOfS <X) CMI 

en *> >H PO u> oj t-< 


-4 

v-minoK 
rwOo-Ho 


HCTC 


ncrcp 
^hoin 


^N-O^CTHr 


loin 

OCM 
HH 


X 

ep 

o 

U f 

> ^ 


"jSSS5S 


t-CMro 




81 

(M 

M 
H 

to 
to 


CTeov 


jfOCT 


ssls?: 


OUD 
-VJO 


f LARVAE THIRO INSTAR LAF 
TP ^n.nny MFW FniTfir.F PHP nnnr.i AQ FTP 


: M 

1- 


> 

;" 

1 


CPoO 
^* 
CM 

rv. jr 


CT^CTecjWCPffJcCfs. 


^00 


CM. 

N. jh< 

JDPOv 
OPOV 


KJN^- 

ainro 

3KiO 
3fOCD 


z 


ro> 

oj- 

^vD 
fK 


^ CP t N. cr cio * cn CP a* 

; 1 
8CLnr"iN.oiroovjDiHvD 


0-HC 

Hin 


HCM 
aeno 


CMCMI- 
*CMC 


, 

voor 

f'O'Hd 
00( 


fa, 

3O 
30 


Hin 


H c\> oo CM ro r*> Jf" In uj 

, 


oOfOf 


^-H-H 
4CM J? 1 

.^f H 

or^(> 


LTJHJ 


jf'CTl 
CMvD< 


NCM 

rcr 
-to 


"*"* 


HCM3-HC\K 

^j-jte^cM^ 


H ; 

i 

fCDN-CTflO 

-HvDiMoO 


HN,. 
OCMl 


^^MJ^ 
ovoro 

3-H- 


MH 

CMHC 
MCMtV 


CM 

ccjj-v 
CsJPO. 


MOO 
DPs* 
tin 


VJDCT 

SS 


HCVJOv 

Lr>H^ 

rj-H 1 *- 


ICMCMC 

CM-HC 
OfO*. 
CM CMC 


FPO so 

JCM 

3S^ 


XJCO 
JDvjD 


CM. 
OCML 


J-C33CM 

H 

hcnro 
rovO 


tf-UDOC 
HH T 

O^CVJlf 
y-HJ 


CM CMI 

CMN.V 
N.CMI 


*-in 
pin 

Oj- 


V 

<U 

r c 

- c 

3 Qt 

Li. U 


; 

; 

to 

<c 


CM 
OCM 


T-|HH, 

tf'CT'^OCOl 
nC9vj CPH, 


<CMCVI 


\JCM 

"OK5 

join 




H-H 


HCMCs, 








-' 


HM I 


_ _ 




HOJf 


"-*" 11 


LOK-OC 


cro. 


-<OJ 
HH 


r^ar 
H-H 


H-H-r 


HHCSJCMCVJ 


[MCM 



2 fOvfjCMyDfOfOHCMHCpOO 



j^u 

s 

*H 

LJ 

P 

5* 



' k^ ! ;>.' 

p ; ! HC 

w r i i 

ii fa i 

1 O. C I 

: - 'a ! j 

y t-f tn^crerlD^irjcr j-rpooo 



|-in-.^ru'H tvi"- 1 UJ-jrcpw' 
O<J -H CM [JMn If) ) ^ cp ^-4 ^ 

I J 1 HH 

! ' i 



UJ 

'o 



IS 



|cn 
tu !ac 



I 



The mission of the PACIFIC NORTHWEST FOREST AND 
RANGE EXPERIMENT STATION is to provide the knowl- 
edge, technology, and alternatives for present and future 
protection, management, and use of forest, range, and related 
environments. 

Within this overall mission, the Station conducts and 
stimulates research to facilitate and to accelerate progress 
toward the following goals: 

1. Providing safe and efficient technology for inventory, 

protection, and use of resources. 

2. Developing and evaluating alternative methods and levels 
of resource management. 

3. Achieving optimum sustained resource productivity 

consistent with maintaining a high quality forest 
environment. 

The area of research encompasses Oregon, Washington, 
Alaska, and, in some cases, California, Hawaii, the Western 
States, and the Nation. Results of the research are made 
available promptly. Project headquarters are at: 



Anchorage, Alaska La Grande, Oregon 

Fairbanks, Alaska Portland, Oregon 

Juneau, Alaska Olympia, Washington 

Bend, Oregon Seattle, Washington 

CorvaSSis, Oregon Wenatchee, Washington 



Mailing address: Pacific Northwest Forest and Range 

Experiment Station 
809 N.E. 6th Ave. 
Portland, Oregon 97232 



CARNEGIE 
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The FOREST SERVICE of the U.S. Department of Agriculture is dedi 
to the principle of multiple u$e management of the Nation's forest resc 
for sustained yields of wood, water, forage, wildlife, and recre 
Through forestry research, cooperation with the States and private 
owners, and management of the National Forests and National Grasslai 
strives as directed by Congress to provide increasingly greater serv 
a growing Nation. 

The U.S. Department of Agriculture Is an Equal Opportunity Emp 

Applicants for all Department programs will be given equal considei 
without regard to age, race, color, sex, religion, or national origin. 



